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Abstract 

The approach to the theory of many-particle interacting systems from a unified standpoint, 
based on the variational principle for free energy is reviewed. A systematic discussion is given 
of the approximate free energies of complex statistical systems. The analysis is centered around 
the variational principle of N. N. Bogoliubov for free energy in the context of its applications to 
various problems of statistical mechanics and condensed matter physics. The review presents a 
terse discussion of selected works carried out over the past few decades on the theory of many- 
particle interacting systems in terms of the variational inequalities. It is the purpose of this 
paper to discuss some of the general principles which form the mathematical background to 
this approach, and to establish a connection of the variational technique with other methods, 
such as the method of the mean (or self-consistent) field in the many-body problem, in which 
the effect of all the other particles on any given particle is approximated by a single averaged 
effect, thus reducing a many-body problem to a single-body problem. The method is illustrated 
by applying it to various systems of many-particle interacting systems, such as Ising and 
Heisenberg models, superconducting and superfluid systems, strongly correlated systems, etc. 
It seems likely that these technical advances in the many-body problem will be useful in 
suggesting new methods for treating and understanding many-particle interacting systems. 
This work proposes a new, general and pedagogical presentation, intended both for those who 
are interested in basic aspects, and for those who are interested in concrete applications. 
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1 Introduction 


The fundamental works of N.N. Bogoliubov on many-body theory and quantum held theory Eca 
OS], on the theory of phase transitions, and on the general theory of interacting systems provided 
a new perspective in various fields of mathematics and physics. The variational principle of N. N. 
Bogoliubov hieieisus] is a useful working tool and has been widely applied to many problems 
of physical interest. It has a well-established place in the many-body theory and condensed 
matter physics [a eubuhi uni nu nil H3i mj. The variational principle of N. N. Bogoliubov has 
led to a better understanding of various physical phenomena such as superfluidity e 0 m m, 
superconductivity E El El n m, phase transitions E ei s si na fl6lj and other cooperative 
phenomena E [ED ED EO , etc. 

Variational methods in physics and applied mathematics were formulated long ago P2E3EDE2 
1231 [Ml 12511261 [271 [25] . It was Maupertuis [25], who wrote in 1774 the celebrated statement: 

Nature, in the production of its effects, does so always by simplest means. 

Since that time variational methods have become an increasingly popular tool in mechanics, hy¬ 
drodynamics, theory of elasticity, etc. Moreover, the variational methods are useful and workable 
tools for many areas of the quantum theory of atoms and molecules EDESJEQIEDEZ], statistical 
many-particle physics and condensed matter physics. The variational methods have been applied 
widely in quantum-mechanical calculations [2T1 132] l29 l l30 l l3 p 132] . in theory of many-particle inter¬ 
acting systems E 0 0 13113 ED (HI 13 M and theory of transport processes [33j El] • As a result 
of these efforts many important and effective methods were elaborated by various researchers. 
Prom the other hand the study of the quasiparticle excitations in many-particle systems has been 
one of the most fascinating subjects for many years EaEunuiH]. The quantum field-theoretical 
techniques have been widely applied to the statistical treatment of a large number of interacting 
particles. Many-body calculations are often done for model systems of statistical mechanics using 
the perturbation expansion. The basic procedure in many-body theory is to find the relevant 
unperturbed Hamiltonian and then take into account the small perturbation operator. This pro¬ 
cedure, which works well for the weakly interacting systems, needs the suitable reformulation for 
the many-body systems with complicated spectra and strong interaction. 

The considerable progress in studying the spectra of elementary excitations and thermodynamic 
properties of many-body systems has been for the most part due to the development of the 
temperature-dependent Green functions methods EaEunuE]. The very important concept of 
the whole method is the concept of the generalized mean field HD E3 E3 ED- These generalized 
mean fields have a complicated structure for the strongly correlated case and are not reduced 
to the functional of the mean densities of the electrons. The concept of the generalized mean 
fields and the relevant algebra of operators from which the corresponding Green functions are 
constructed are the central ones to our treatment of the strongly interacting many-body systems. 
It is the purpose of this paper to discuss some of the general principles which form the physical 
and mathematical background to the variational approach, and to establish the connection of the 
variational technique with other methods in the theory of many-body problem. 

2 The Variational Principles of Quantum Theory 

It is well known that in quantum mechanics the eigenfunction ipi of the lowest state of any system 
has the property of making the integral 

J ip* HipiCpr (1) 
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a minimum. The value of integral is the corresponding eigenvalue Ej of the Hamiltonian H of a 
system. These circumstances lead to a specific approximate method (the variational method) of 
finding ipi and £) by minimizing integral JT]) among a restricted class of functions. 

The variational method [2T> ;32) 129], 130] [31] enables one to make estimates of energy levels by using 
trial wave functions ipT 

= f^Hip T d 3 r 

T f i/j^Td 3 r 


The ground state Eq gives the lowest possible energy the system can have. Hence, for the approx¬ 
imation of the ground state energy one would like to minimize the expectation value of the energy 
with respect to a trial wave function. 

In other words the variational principle states that the ground-state energy of a quantum mechan¬ 
ical system is less than or equal to the expectation value of the Hamiltonian with an arbitrary 
wave function. Given a trial wave function with adjustable parameters, the best values of the 
parameters are those which minimize the expectation value of the Hamiltonian. The variational 
principle consists in adjusting the available parameters, so as to maximize this lower bound. 

An important method of finding approximate ground state energies and wave functions is called 
by the Rayleigh-Ritz variational principle [21 , 1291 [30l [3T] . The Rayleigh-Ritz variational principle 
for the ground state energy is the starting point of many computations and approximations in 
quantum mechanics and quantum chemistry of atoms and molecules. This principle states that 
the expectation value of H in any state | ijj) is always greater than or equal to the ground state 
energy, E 0 


(mm 

(V#) 


> Eo, 


(3) 


or 


(H) > MHW) > E 0 . 


(4) 


Here \ijj) € Q is arbitrary pure quantum state and H is a Hamiltonian acting on a Hilbert space 
Q. This relation becomes equality only when i/j = i/jq. Thus this principle gives the upper bound 
to the ground state energy. 

It will be instructive also to remind how the variational principle of quantum mechanics com¬ 
plements the perturbation theory |38], [39|. For this aim let us consider the Rayleigh-Schrodinger 
perturbation expansion. The second-order level-shift E® of the ground state of a system have the 
form 


*§ = £ 
3 +0 


(r\v\^)(^\v\r 

(EP-Ej) 


E 

3 +0 


I Hr 


Oil 


(E*-E 3 y 


(5) 


where Voj = Voj\4>j) and |^°) is the unperturbed ground state. It is clear then that E® is 
always negative. 

The variational principle of quantum mechanics states that the ground state energy E° for the 
total Hamiltonian H is the minimum of the energy functional 


e {*} = (mm, ( 6 ) 

where T is a trial wave function. It should be noted that it is possible to establish that the sum 
of all the higher-order level shifts E®, starting with n = 2, will be negative, providing that the 
relevant perturbation series will converges to E°. 

To confirm this statement let us consider again the Hamiltonian 

H = H 0 + XV. (7) 
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It is reasonable to suppose that the ground state energy E° = E°( A) and the ground state 
T = W(A) of the Hamiltonian H are analytic functions (at least for small A). Note that when 
one considers the many-body problem, the concept of relative boundedness is of use, where a 
perturbation XV is small compared to H o in the sense that 

(H 0 ) 2 > (X 2 V 2 ). 

This means simply that the eigenvalues of the operator ((Ho) 2 — (X 2 V 2 )) are non-negative. 

Then the corresponding perturbation expansion may be written in the form 

Eo = 4 0) + A^ 1} + A 2 4 2) + a3 4 3) + • • •, (8) 

where Eq = (fj°\H\il> 0 ) and Eq = ('i/>°|R| 

The variational approach states that 

E 0 = mm((^\H 0 + XV\^)). (9) 

Thus we obtain 

A 2 E {2) + A 3 4 3) + ... = E 0 - (4 0) + a4 1} ) = 

(mm{(T>\H 0 + AH|T)} - </| H 0 + AR|V>°>) . (10) 

In this expression the second part must satisfy the condition 

(mm{(^\H 0 + AR|T)} - {ip°\H 0 + AH|^ 0 )) < 0. (11) 

In addition, in general case the relevant ground state T which yields a minimum will not coincide 
with if 0 . Thus we obtain 

X 2 E {2) + A 3 £j 3) + ... < 0. (12) 

The last inequality can be rewritten as 

E {2) < (A+ A 2 4 4) + • • •)■ (13) 

( 2 ) 

In the limit A —>• 0 we have that Eq < 0. Thus the variational principle of quantum mechanics 
confirms the results of the perturbation theory [3D] . 

It is worth mentioning that the Rayleigh-Ritz variational method has a long and interesting his¬ 
tory BH S3 S3]. Rayleigh’s classical book Theory of Sound was first published in 1877. In it 
are many examples of calculating fundamental natural frequencies of free vibration of continuum 
systems (strings, bars, beams, membranes, plates) by assuming the mode shape, and setting the 
maximum values of potential and kinetic energy in a cycle of motion equal to each other. This 
procedure is well known as Rayleigh’s Method. In 1908, Ritz laid out his famous method for de¬ 
termining frequencies and mode shapes, choosing multiple admissible displacement functions, and 
minimizing a functional involving both potential and kinetic energies. He then demonstrated it in 
detail in 1909 for the completely free square plate. In 1911, Rayleigh wrote a paper congratulating 
Ritz on his work, but stating that he himself had used Ritz’s method in many places in his book 
and in another publication. 

Subsequently, hundreds of research articles and many books have appeared which use the method, 
some calling it the ’’Ritz method” and others the ”Rayleigh-Ritz method.” The article [41] exam¬ 
ined the method in detail, as Ritz presented it, and as Rayleigh claimed to have used it. A. W. 
Leissa [H] concluded that, although Rayleigh did solve a few problems which involved minimiza¬ 
tion of a frequency, these solutions were not by the straightforward, direct method presented by 
Ritz and used subsequently by others. Therefore, Rayleigh’s name should not be attached to the 
method. Additional informative comments were carried out in Refs. [42[ 133] 
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3 The Helmholtz Free Energy and Statistical Thermodynamics 

Variational methods in thermodynamics and statistical mechanics have been used widely since J. 
W. Gibbs groundbreaking works [26l [271 [28] . According to Gibbs approach a workable procedure 
for the development of the statistical mechanical ensemble theory is to introduce the Gibbs entropy 
postulate. Hence, as a result of the Gibbs ensemble method, the entropy S can be expressed in 
the form of an average for all the ensembles, namely, 

S(N,V,E ) = -kB^Piln-Pi = = k B lnfi(JV, V, E), (14) 

l 

where the summation over i denotes a general summation over all states of the system and pi is 
the probability of observing state i in the given ensemble and hs is the Boltzmann constant. This 
relation links entropy S and probability jy. 

It can be said that in this context the entropy is a state function which is according to the second 
law mm is defined by the relation 


dS = (3(dE - dF). (15) 

The energy E and the Helmholtz free energy F are the state functions [23 04]. The proportionality 
coefficient f5 was termed as the thermodynamic temperature (/3 = 1 /ksT) of the surrounding with 
which the system exchanges by heat Q and work W. 

Thus the postulate of equal probabilities in the microcanonical ensemble [45] and the Gibbs 
entropy postulate can be considered as a convenient starting points for the development of the 
statistical mechanical ensemble theory in a standard approach m- 

After postulating the entropy by means of Eq. (11411 , the thermodynamic equilibrium ensembles are 
determined by the following criterion for equilibrium: 

(8S)e,v,n = 0- (16) 

This variational scheme is used for each ensemble (microcanonical, canonical and grand canonical) 
with different constraints for each ensemble. In addition, this procedure introduces Lagrange 
multipliers which, in turn, must be identified with thermodynamic intensive variables (T, P). From 
the other hand, the procedure of introducing Lagrange multipliers and the task of identifying them 
with the thermodynamic intensive properties can be clarified by invoking a more general criterion 
for thermodynamic equilibrium. 

From the Gibbs entropy postulate, Eq. (ll4D . the definitions of average and the normalization 
constraint J2iPi = 1 one obtains 


5S = + \npi)5pi, (17) 

i 

5E = Y j Eidpi , (18) 

i 

6V = Y, Vi6 p^ ( 1£ ») 

i 

^2 6 Pi = °- ( 2 °) 

i 

Using a Lagrange multiplier A together with the variational condition, we obtain 

'y \ {Ej A EVi + A + ksT + kBT\npi)5pi > 0. (21) 

i 
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Here all 5pi are considered as the independent variables. Thus we deduce that 

p i = exp(-p\-l-P(PV i + E i )'), /3 = (k B T)-\ (22) 

The Lagrange multiplier A, can be determined directly from the definition of entropy (1141) . 


S = —k B Pi ( 


/Ei + PVi + A + k B T\ _ (E + PV + A + k B T^j 
V kPT ) ~ T 


(23) 


Thus we arrive at 


A + k B T = TS - E - PV = -G, (24) 

Pi = expP(G-PVi-El). (25) 

Here G is the Gibbs energy (or Gibbs free energy). It may also be defined with the aid of the 
Helmholtz free energy G = H — TS. Here H(S, P, N ) is the enthalpy [ 341 . 

The usefulness of the thermodynamic potentials G and F may be clarified within the statistical 
thermodynamics m For the microcanonical ensemble one should substitute Ei = E and V* = V, 
which are fixed for every system and since G — PV — E = S, Eo. (l25l) becomes 

Vi = e ~ s / kB . (26) 


For the canonical ensemble one should substitute Vi = V, which is given for each system and in 
this case Eg. (1251) can be written as 


Pi = e^ F ~ Ei \ 


(27) 


Here F = G — PV denotes the Helmholtz free energy. Thus the free energy F is defined by 

F = E-TS. (28) 


The Helmholtz free energy describes an energy which is available in the form of useful work. 

The second law of thermodynamics asserts that in every neighborhood of any state A in an 
adiabatically isolated system there exist other states that are inaccessible from A. This statement 
in terms of the entropy S and heat Q can be formulated as 


dS = dQ/T + da. (29) 

Thus the only states available in an adiabatic process ( dQ = 0, or dS = da) are those which 
lead to an increase of the entropy S. Here da > 0 defines the entropy production a due to the 
irreversibility of the transformation. 

It is of use to analyze the expression 

dF = dE — TdS - SdT = -SdT - Tda - PdV + PiN t . (30) 

Free energy change A F of a system during a transformation of a system describes a balance of the 
work exchanged with the surroundings. If A F > o, A F represents the minimum work that must 
be incurred for the system to carry out the transformation. In the case A F < o, | AT 1 1 represents 
the maximum work that can be obtained from the system during the transformation. It is obvious 
that 

dF = dE- Tda - SdT. (31) 


6 




In a closed system without chemical reaction and in the absence of any other energy exchange, 
the variation A F = —SdT — TdS — PdV + /M,Ni can be rewritten in the form 

dF = -Tda < 0. (32) 

It means that function F decreases and tends towards a minimum corresponding to equilibrium. 
Thus the Helmholtz free energy is the thermodynamic potential of a system subjected to the 
constraints constant T,V,Ni. 

The Gibbs free energy (free enthalpy) is defined by 


G = H - TS = F + PV. (33) 

The physical meaning of the Gibbs free energy is clarified when considering evolution of a system 
from a certain initial state to a final state. The Gibbs free energy change AG then represents 
the work exchanged by the system with its environment and the work of the pressure forces, 
during a reversible transformation of the system. Here H = E + VP = TS + VP + )G /ijiVj is 
the thermodynamic potential of a system termed by enthalpy [33] , The Gibbs free energy is the 
thermodynamic potential of a system subjected to the constraints constant T, P, N^. In this case 

dG = -Tda < 0. (34) 

Thus in the closed system without chemical reaction and in the absence of any other energy 
exchange at constant temperature, pressure and amount of substance the function G can only 
decrease and reach a minimum at equilibrium. 

It will be of use to mention another class of thermodynamic potentials, termed by the Massieu- 
Planck functions. These objects may be deduced from the fundamental relations in the entropy 
representations, S = S(E,V,N). The corresponding differential form may be written as 

dS = ±dE + ^dV-^dN. (35) 

Thus the suitable variables for a Legendre transform will be 1/T, P/T and p/T. In some cases 
working with these variables is more convenient. 

It is worth noting that in terms of the Gibbs ensemble method the free energy is the thermo¬ 
dynamic potential of a system subjected to the constraints constant T. V. Ni. Moreover, the 
thermodynamic potentials should be defined properly in the thermodynamic limit. The problem 
of the thermodynamic limit in statistical physics was discussed in detail by A. L. Kuzemsky [45]. 
To clarify this notion, let us consider the logarithm of the partition function Q(0, V, N) 

F(8,V,N) = —Q\nQ(0,V,N). (36) 

This expression determines the free energy F of the system on the basis of canonical distribu¬ 
tion. The standard way of reasoning in the equilibrium statistical mechanics do not requires the 
knowledge of the exact value of the function F(8,V,N). For real system it is sufficient to know 
the thermodynamic (infinite volume) limit [15 1 [271 - 351 [35] 

Jim F ( e '£ N ) \ = m V/N). (37) 

JV-KX) A ' 

Here f(8,V/N ) is the free energy per particle. It is clear that this function determines all the 
thermodynamic properties of the system. 

Let us summarize the criteria for equilibrium briefly. In a system of constant V and S, the 
internal energy has its minimum value, whereas in a system of constant E and V, the entropy has 
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its maximum value. 

It should be noted that the pair of independent variables (V, S) is not suitable one because of the 
entropy is not convenient to measure or control. Hence it would be of use to have fundamental 
equations with independent variables that is easier to control. The two convenient choices are 
possible. First we take the P and T pair. From the practical point of view this is a convenient 
pair of variables which are easy to control (measure). For systems with constant pressure the best 
suited state function is the Gibbs free energy (also called free enthalpy) 

G = H — TS. (38) 

Second relevant pair is V and T. For systems with constant volume (and variable pressure), the 
suitable suited state function is the Helmholtz free energy 

F = E — TS. (39) 

Any state function can be used to describe any system (at equilibrium, of course), but for a given 
system some are more convenient than others. The change of the Helmholtz free energy can be 
written as 

dF = dE — TdS - SdT. (40) 

Combining this equation with dU = TdS — PdV we obtain the relation of the form 


dF = -PdV - SdT. 


In terms of variables (T, V) we find 


: dF 

dF=t er 


„ /dF 
dT + ( — 


Comparing the equations one can see that 

'dl 

.dT 


dF 

S = -(^= 


v \dV 


P= (dF 

v’ \dV 


dV. 


(41) 


(42) 


(43) 


At constant T and V the equilibrium states corresponds to the minimum of Helmholtz free energy 
(dF = 0). From F = E — TS we may suppose that low values of F are obtained with low values 
of E and high values of S. 

In terms of a general statistical-mechanical formalism m m m, a many-particle system with 
Hamiltonian H in contact with a heat bath at temperature T in a state described by the statistical 
operator p has a free energy 

F = Tr(pH)+k B TTr(plnp). (44) 

The free energy takes its minimum value 


F eq = -k B T In Z (45) 

in the equilibrium state characterized by the canonical distribution 

p eq = Z~ x exp(-Hf3)- Z = Tr exp(—if/3). (46) 

Before turning to the next topic, an important remark about the free energy will not be out of place 
here. I. Novak m attempted to give a microscopic description of Le Chatelier’s principle [38| in 
statistical systems. Novak have carried out interesting analysis based on microscopic descriptors 
(energy levels and their populations) that provides visualization of free energies and conceptual 
rationalization of Le Chatelier’s principle. The misconception ” nature favors equilibrium ” was 
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highlighted. This problem is a delicate one and requires a careful discussion |49j. Dasmeh et al. 
showed [49] that Le Chatelier’s principle states that when a system is disturbed, it will shift its 
equilibrium to counteract the disturbance. However for a chemical reaction in a small, confined 
system, the probability of observing it proceed in the opposite direction to that predicted by 
Le Chatelier’s principle, can be significant. Their study provided a molecular level proof of Le 
Chatelier’s principle for the case of a temperature change. Moreover, a new, exact mathematical 
expression was derived that is valid for arbitrary system sizes and gives the relative probability 
that a single experiment will proceed in the endothermic or exothermic direction, in terms of a 
microscopic phase function. They showed that the average of the time integral of this function 
is the maximum possible value of the purely irreversible entropy production for the thermal 
relaxation process. The results obtained were tested against computer simulations of the unfolding 
of a polypeptide. It was proven that any equilibrium reaction mixture on average responds to a 
temperature increase by shifting its point of equilibrium in the endothermic direction. 


4 Approximate Calculations of Helmholtz Free Energy 

Statistical mechanics provides effective and workable tools for describing the behavior of the sys¬ 
tems of many interacting particles. One of such approaches for describing systems in equilibrium 
consists in evaluation the partition function Z and then the free energy. 

Now we must take note of the different methods for obtaining the approximate Helmholtz free en¬ 
ergy in the theory of many-particle systems. Roughly speaking, there are two approaches, namely 
the perturbation method and the variational method. 

Thermodynamic perturbation theory [501 EH ISZl [53] may be applied to systems that undergo a 
phase transition. It was shown [54] that certain conditions are necessary in order that the applica¬ 
tion of the perturbation does not change the qualitative features of the phase transition. Usually, 
the shift in the critical temperature are determined to two orders in the perturbation parameter. 
Let us consider here the perturbation method [54] very briefly. 

In paper [54] authors considered a system with Hamiltonian Hq , that undergoes a phase transition 
at critical temperature T®,. The task was to determine for what class of perturbing potentials V 
will the system with Hamiltonian Hq + V have a phase transition with qualitatively the same 
features as the unperturbed system. In their paper authors [54J studied that question using ther¬ 
modynamic perturbation theory [50] ED [52]. They found that an expansion for the perturbed 
thermodynamic functions can be term by term divergent at the critical temperature Tq for a class 
of potentials V. Under certain conditions the series can be resummed, in which case the phase 
transition remains qualitatively the same as in the unperturbed system but the location of the 
critical temperature is shifted. 

The starting point was the partition function Zq for a system whose Hamiltonian is Hq 


Zq = Tr exp(—LT 0 /3). 


(47) 


For a system with Hamiltonian Hq + AU, the partition function Z is given by 


Z = Trexp[—(iL 0 + XV)/3\. 


(48) 


It is possible to obtain formally an expansion for Z in terms of the properties of the unperturbed 
system by expanding that part of the exponential containing the perturbation in the following 
way [54] when V and Hq commute: 

OO OO 1 

z = Tr(exp(—tfo/3) £ ^[-A P] n V n ) = Z 0 £ -[-A /3} n (V n ) 0 , (49) 

n n 
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where 


(50) 


Z o (V n ) o = Ti(e X .p[-H o 0\V n ). 

Then the expression for Z can be written as 

7 00 1 

-=exp(m + £-(-Affl"<F”}„]) (51) 

^ n 

The free energy per particle / is given by 

1 00 

Pf, = f>h - Jf ln(l + £ ^(-W( V) 0 ) • (52) 

n 

where f p and /o are the perturbed and unperturbed free energy per particle, respectively, and N 
is the number of particles in the system. The standard way to proceed consists of in expanding 
the logarithm in powers of A. As a result one obtains [53] 

Pfp = Pfo + jjr(V) 0 - ^^(V>0 - (V)l) 
d ((^ 3 )° — 3(R 2 )o(R)o + (^)o) + • • • (53) 

To proceed, it is supposed usually that the thermodynamics of the unperturbed system is known 
and the perturbation series (if they converge) may provide us with suitable corrections. If the 
terms in the expansion diverge, they may, in principle, be regularized under some conditions. For 
example, perturbation expansions for the equation of state of a fluid whose intermolecular poten¬ 
tial can be regarded as consisting of the sum of a strong and weak part give reasonable qualitative 
results [53 [56]. 

In paper by Fernandes [53] author investigated the application of perturbation theory to the 
canonical partition function of statistical mechanics. The Schwinger and Rayleigh-Schrodinger 
perturbation theory were outlined and plausible arguments were formulated that both should give 
the same result. It was shown that by introducing adjustable parameters in the unperturbed or 
reference Hamiltonian operator, one can improve the rate of convergence of Schwinger perturba¬ 
tion theory. The same parameters are also suitable for Rayleigh-Schrodinger perturbation theory. 
Author discussed also a possibility of variational improvements of perturbation theory and gave 
a simpler proof of a previously derived result about the choice of the energy shift parameter. It 
was also shown that some variational parameters correct the anomalous behavior of the partition 
function at high temperatures in both Schwinger and Rayleigh-Schrodinger perturbation theories. 
It should be stressed, however, that the perturbation method is valid for small perturbations only. 
The variational method is more flexible tool [381 [391 EZl ESI |59l EOl ED [62] and in many cases is 
more appropriate in spite of the obvious shortcomings. But both the methods are interrelated 
deeply [39] and enrich each other. 

R. Peierls mmm pointed at the circumstance that for a many-particle system in thermal 
equilibrium there is a minimum property of the free energy which may be considered as a gen¬ 
eralization of the variational principle for the lowest eigenvalue in quantum mechanics. Peierls 
attracted attention to the fact that the free energy has a specific property which can be formu¬ 
lated in the following way. Let us consider an arbitrary set of orthogonal and normalized functions 
{<^i, <£> 2 , • • • (^n, • • The expectation value of the Hamiltonian H for nth of them will be written 
as 

Hnn = J VnHVndr. (54) 
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The statement is that for any temperature T the function 


F = -k B T log Z = -k B T log ^ exp [-H nn /3\ 


( 55 ) 


which would represent the free energy if the H nn were the true eigenvalues, is higher than the 
true free energy 

F 0 = -k B T log Z 0 = —k B T log exp[—E n (3\, (56) 


or 


F > F 0 . (57) 

This is equivalent to saying that the partition function, as formed by means of the expectation 
values H nn : 

Z = exp [-H nn j3] (58) 


is less than the true partition function 


Z 0 = y^exp [~E n /3\, 


(59) 


or 


Z 0 = y~^exp [~E n /3] > Z = y~^exp [-H nn /3\. 


(60) 


Peierls [63j formulated the more general statement, namely, that if f(E ) is a function with the 
properties 

df n d 2 f 


the expression 


is less than 


dE < °’ dE 2 > °’ 


f = Y J f{H nn ) 

n 

f 0 = J2f(E n ). 


(61) 

(62) 

(63) 


To summarize, Peierls has proved a kind of theorem a special case of which gives a lower bound to 
the partition sum and hence an upper bound to the free energy of a quantum mechanical system 


^exp [~E k /3] > y^exp [-H nn /3]. 


(64) 


When j3 —>• oo the theorem is obvious, reducing to the fundamental inequality Ek < H nn for all 
n. However for finite (3 it is not so obvious since higher eigenvalues of H do not necessarily lie 
lower than corresponding diagonal matrix elements H nn . T. D. Schultz [64] skilfully remarked 
that, in fact, the Peierls inequality does not depend on the fact that exp[— E/3\ is a monotonically 
decreasing function of E, as might be concluded from the original proof. It depends only on the 
fact that the exponential function is concave upward. Schultz [6l] proposed a simple proof of the 
theorem under this somewhat general condition. 

Let (p n be a complete orthonormal set of state vectors and let A be an Hermitean operator which 
for convenience is assumed to have a pure point spectrum with eigenvalues ak and eigenstates 
Let f(x ) be a real valued function such that 


dx 2 


> 0 


(65) 
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in an interval including the whole spectrum of a k - Then if Tr/(A) exists it can be proven that 
the following statement holds 

Tr/(A) > ^ f( a nn), (66) 

n 

where a nn = (n|A|n). The equality holds if and only if the <p n are the eigenstates of A. 

Since 

Tr/(A) = J>|/(A)|n) (67) 

n 

it is sufficient for the proof to point out that the relation (1661) follows from 

(n|/(A)|n) > f(a nn ), (68) 

which is valid for all n. The inequalities (1681) were derived from 

f(a k ) > /( ^nn ) + (ofc ~ ^nn)/ (^nn ), (69) 

which is a consequence of equation (|65l) . the right hand side for fixed n being the line tangent 
to /(a/c) at a nn . Multiplying (|69l) by |(n|fc)| 2 and summing on k one obtains (f68l) . Schultz j6l j 
observed further that the equality in (f68l) holds if and only if |(n|fc)| 2 = 0 unless a k = a nn , i.e. if 
and only if ip n is an eigenstate of A. 

If /(A) is positive definite, then the set (p n need not be complete, since the theorem is true even 
more strongly if positive terms are omitted from the sum f(a nn ). 

With the choice /(A) = exp(—A) and A = H/3, the original theorem of Peierls giving an upper 
bound to the free energy is reproduced. With A = (H — ^N)/3 we have an analogous theorem for 
the grand potential. The theorem proved by Schultz [64] is a generalization in that it no longer 
requires f(x ) to be monotonic; it requires only that Tr/(A) be finite which can occur even if /(x) 
is not monotonic provided A is bounded. 

Peierls variational theorem was discussed and applied in a number of papers (see, e.g. Refs. [6411651 
[66l [67] ). It has much more generality than, say, the A. Lidiard [68] consideration on a minimum 
property of the free energy. Lidiard [68] derived the approximate free energy expression in a way 
which shows a strong analogy with the approximate Hartree method of quantum mechanics. By 
his derivation he refined the earlier calculations made by Koppe an Wohlfarth in the context of 
description of the influence of the exchange energy on the thermal properties of free electrons in 
metals. 

5 The Mean Field Concept 

In general case, a many-particle system with interactions is very difficult to solve exactly, except 
for special simple cases. Theory of molecular (or mean) field permits one to obtain an approximate 
solution to the problem. In condensed matter physics, mean field theory (or self-consistent field 
theory) studies the behavior of large many-particle systems by studying a simpler models. The 
effect of all the other particles on any given particle is approximated by a single averaged effect, 
thus reducing a many-body problem to a single-body problem. 

It is well known that molecular fields in various variants appear in the simplified analysis of many 
different kinds of many-particle interacting systems. The mean field concept was originally for¬ 
mulated for many-particle systems (in an implicit form) in Van der Waals [69, 70] dissertation 
”On the Continuity of Gaseous and Liquid States”. Van der Waals conjectured that the volume 
correction to the equation of state would lead only to a trivial reduction of the available space 
for the molecular motion by an amount b equal to the overall volume of the molecules. In reality, 


12 


the measurements lead him to be much more complicated dependence. He found that both the 
corrections should be taken into account. Those were the volume correction b, and the pressure 
correction a/V 2 , which led him to the Van der Waals equation [70] . Thus, Van der Waals came 
to conclusion that ’’the range of attractive forces contains many neighboring molecules”. The 
equation derived by Van der Waals was similar to the ideal gas equation except that the pressure 
is increased and the volume decreased from the ideal gas values. Hence, the many-particle be¬ 
havior was reduced to effective (or renormalized) behavior of a single particle in a medium (or a 
field). The later development of this line of reasoning led to the fruitful concept, that it may be 
reasonable to describe approximately the complex many-particle behavior of gases, liquids, and 
solids in terms of a single particle moving in an average (or effective) field created by all the other 
particles, considered as some homogeneous (or inhomogeneous) environment. 

Later, these ideas were extended to the physics of magnetic phenomena mmmm, where 
magnetic substances were considered as some kind of a specific liquid. This approach was elab¬ 
orated in the physics of magnetism by P. Curie and P. Weiss. The mean field (molecular field) 
replaces the interaction of all the other particles to an arbitrary particle [73] . In the mean field 
approximation, the energy of a system is replaced by the sum of identical single particle energies 
that describe the interactions of each particle with an effective mean field. 

Beginning from 1907 the Weiss molecular-field approximation became widespread in the theory of 
magnetic phenomena nzi Ei eh [72], and even at the present time it is still being used efficiently. 
Nevertheless, back in 1965 it was noticed that 131 = 

” The Weiss molecular field, theory plays an enigmatic role in the statistical mechanics 
of magnetism ”. 

In order to explain the concept of the molecular field on the example of the Heisenberg ferromagnet 
one has to transform the original many-particle Hamiltonian 

H = j) SiSj - gfiB H Y s i- ( 7 °) 

ij i 

into the following reduced one-particle Hamiltonian 

n = - 2 p 0 p B s- h {mf) . 

The coupling coefficient J{i — j) is the measure of the exchange interaction between spins at the 
lattice sites i and j and is defined usually to have the property J(i - j = 0) = 0. This transformation 
was achieved with the help of the identity 

S ■ S' = S ■ {S') + (5) • S' - (5) • {S') + C. 

Here, the constant C = (S — (S)) ■ ( S' — {S')) describes spin correlations. The usual molecular- 
field approximation is equivalent to discarding the third term in the right hand side of the above 
equation, and using the approximation C ~ {C) = {S ■ S') — {S) ■ {S'), for the constant C. 

There is large diversity of the mean-field theories adapted to various concrete applications mu 

EH [72]. 

Mean field theory has been applied to a number of models of physical systems so as to study 
phenomena such as phase transitions |75lE6|. One of the first application was Ising model mm 
EDE2]. Consider the Ising model on an IV-dimensional cubic lattice. The Hamiltonian is given 

by __ 

H = -jYSiSj-hYSi, (71) 

(*j) * 
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where the Yluj) indicates summation over the pair of nearest neighbors (i. j) , and S t = ±1 and 
Sj are neighboring Ising spins. A. Bunde m has shown that in the correctly performed molecular 
held approximation for ferro- and antiferromagnets the correlation function (S'(q)S'(-q)) should 
fulfill the sum rule 

iV-l]T<S(q)<?(-q)) = l. (72) 

<? 

The Ising model of the ferromagnet was considered EZ3 and the correlation function (5(q)5(—q)) 
was calculated 


<S(q)S(-q)>= N~ l ]T 


1 


-1 


1 


1 — /3«7(q) ’ 


(73) 


which obviously fulhlls the above sum rule. The Ising model and the Heisenberg model were two 
the most explored models for the applications of the mean held theory. 

It is of instruction to mention that earlier molecular-held concepts described the mean-held in 
terms of some functional of the average density of particles (n) (or, using the magnetic terminology, 
the average magnetization (M)), that is, as F[(n), (M) ]. Using the modern language, one can say 
that the interaction between the atomic spins Si and their neighbors can be equivalently described 
by effective (or mean) held . As a result one can write down 


M.-XolkU’+Z-fU 

The mean held can be represented in the form (in the case T > Tq) 

hW) = Y J j( R j i){s l ). 


(74) 


(75) 


Here, h ext is the external magnetic held, xo is the system’s response function, and J{Rji ) is the 
interaction between the spins. In other words, in the mean-held approximation a many-particle 
system is reduced to the situation, where the magnetic moment at any site aligns either parallel 
or anti-parallel to the overall magnetic held, which is the sum of the applied external held and 
the molecular held. 

Note that only the ” averaged ” interaction with i neighboring sites was taken into account, while 
the fluctuation effects were ignored. We see that the mean-held approximation provides only a 
rough description of the real situation and overestimates the interaction between particles. At¬ 
tempts to improve the homogeneous mean-held approximation were undertaken along different 

directions mummmmm- 

An extremely successful and quite nontrivial approach was developed by L. Neel mumm, 
who essentially formulated the concept of local mean fields (1932). Neel assumed that the sign 
of the mean-held could be both positive and negative. Moreover, he showed that below some crit¬ 
ical temperature (the Neel temperature) the energetically most favorable arrangement of atomic 
magnetic moments is such, that there is an equal number of magnetic moments aligned against 
each other. This novel magnetic structure became known as the antiferromagnetism [5, T7\. 

It was established that the antiferromagnetic interaction tends to align neighboring spins against 
each other. In the one-dimensional case this corresponds to an alternating structure, where an 
”up” spin is followed by a ’’down” spin, and vice versa. Later it was conjectured that the state 
made up from two inserted into each other sublattices is the ground state of the system (in the 
classical sense of this term). Moreover, the mean-field sign there alternates in the ’’chessboard” 
(staggered) order. 

The question of the true antiferromagnetic ground state is not completely clarified up to the 
present time [l7l[5l[7lU72l[35U36l[37]. This is related to the fact that, in contrast to ferromag- 
nets, which have a unique ground state, antiferromagnets can have several different optimal states 
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with the lowest energy. The Neel ground state is understood as a possible form of the system’s 
wave function, describing the antiferromagnetic ordering of all spins. Strictly speaking, the ground 
state is the thermodynamically equilibrium state of the system at zero temperature. Whether the 
Neel state is the ground state in this strict sense or not, is still unknown. It is clear though, 
that in the general case, the Neel state is not an eigenstate of the Heisenberg antiferromagnet’s 
Hamiltonian. On the contrary, similar to any other possible quantum state, it is only some linear 
combination of the Hamiltonian eigenstates. Therefore, the main problem requiring a rigorous 
investigation is the question of Neel state stability HZI- In some sense, only for infinitely large 
lattices, the Neel state becomes the eigenstate of the Hamiltonian and the ground state of the sys¬ 
tem. Nevertheless, the sublattice structure is observed in experiments on neutron scattering 023, 
and, despite certain worries, the actual existence of sublattices is beyond doubt. 

Once Neel’s investigations were published, the effective mean-field concept began to develop at 
a much faster pace. An important generalization and development of this concept was proposed 
in 1936 by L. Onsager [75] in the context of the polar liquid theory. This approach is now called 
the Onsager reaction field approximation. It became widely known, in particular, in the physics 
of magnetic phenomena [7911801 18T] . In 1954, Kinoshita and Narnbu [82J developed a systematic 
method for description of many-particle systems in the framework of an approach which corre¬ 
sponds to the generalized mean-field concept. N.D. Mermin [83] has analyzed the thermal 
Hartree-Fock approximation [S3] of Green function theory giving the free energy of a system not 
at zero temperature. 

Kubo and Suzuki [85j studied the applicability of the mean field approximation and showed that 
the ordinary mean-field theory is restricted only to the region ksT > zJ, where J denotes the 
strength of typical interactions of the relevant system and z the number of nearest neighbours. 
Suzuki [86] has proposed a new type of fluctuating mean-field theory. In that approach the true 
critical point Tq differs from the mean-field value and the singularities of response functions are, 
in general, different from those of the Weiss mean-field theory mm 

Lei Zhou and Ruibao Tao m developed a complete Hartree-Fock mean-field method to study 
ferromagnetic systems at finite temperatures. With the help of the complete Bose transforma¬ 
tion, they renormalized all the high-order interactions including both the dynamic and the kinetic 
ones based on an independent Bose representation, and obtained a set of compact self-consistent 
equations. Using their method, the spontaneous magnetization of an Ising model on a square 
lattice was investigated. The result is reasonably close to the exact one. Finally, they discussed 
the temperature dependences of the coercivities for magnetic systems and showed the hysteresis 
loops at different temperatures. 

Later, various schemes of ’’effective mean-field theory taking into account correlations” were pro¬ 
posed (see Refs. m EZ]). We will see below that various mean-field approximations can be in 
principle described in the framework of the variation principle in terms of the Bogoliubov inequal¬ 
ity e 0 u nni ns] : 


-/r 1 ln(Tre“ /3Hmod ) 


F = —/T 1 \n(Tre-P H ) < 
Xre-^ mod (iL-iL mod ) 

^mod 


(76) 


Here, F is the free energy of the system under consideration, whose calculation is extremely 
involved in the general case. The quantity H mQ d is some trial Hamiltonian describing the effective- 
field approximation. The inequality (1761) yields an upper bound for the free energy of a many- 
particle system. 

It is well known, that the study of Hamiltonians describing strongly-correlated systems is an 
exceptionally difficult many-particle problem, which requires applications of various mathematical 
methods [171 IB5118911901 [91] . In fact, with the exception of a few particular cases, even the ground 
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state of the Hubbard model is still unknown. Calculation of the corresponding quasiparticle 
spectra in the case of strong inter-electron correlations and correct definition of the mean fields 
also turned out to be quite a complicated problem. 

The Hamiltonian of the Hubbard model m is given by: 

H = y^ tija\ a a ja + U/2 ^ n i(7 ni- a . (77) 

ijcr it7 

The above Hamiltonian includes the repulsion of the single-site intra-atomic Coulomb U, and tij , 
the one-electron hopping energy describing jumps from a j site to an i site. As a consequence of 
correlations electrons tend to ’’avoid one another”. Their states are best modeled by atom-like 
Wannier wave functions [</>(r — Rj)\ . The Hubbard model’s Hamiltonian can be characterized by 
two main parameters: U, and the effective band width of tightly bound electrons 

a^-EM 2 ) 172 - 

ij 

The band energy of Bloch electrons e(k) is given by 

e(k ) = N~ 1 Uj exp [—ik(Ri — Rj\, 

k 


where N is the total number of lattice sites. Variations of the parameter 7 = A/f7 allow one to 
study two interesting limiting cases, the band regime (7 1) and the atomic regime (7 —?• 0). 

There are many different approaches to construction of generalized mean-field approximations; 
however, all of them have a special-case character. The method of irreducible Green func¬ 
tions PH (351 USM] allows one to tackle this problem in a more systematic fashion. 

The efficiency of the method of the irreducible Green functions for description of normal and 
superconducting properties of systems with a strong interaction and complicated character of the 
electron spectrum was demonstrated in the papers H3 EH [SEE]. Let us consider the Hubbard 
model m • The properties of this Hamiltonian are determined by the relationship between the 
two parameters: the effective band’s width A and the electron’s repulsion energy U. Drastic 
transformations of the metal-dielectric phase transition’s type take place in the system as the 
ratio of these parameters changes. Note that, simultaneously, the character of the system descrip¬ 
tion must change as well, that is, we always have to describe our system by the set of relevant 
variables. In the case of weak correlation 113 [351 (361 [37] the corresponding set of relevant variables 
contains the ordinary second-quantized Fermi operators and a\ a di a , as well as the number of 
particles operator rii a = a\ a ai a . In the case of strong correlation na [35i iMi [37] the problem is 
highly complicated. 

The Green function in the generalized mean-field’s approximation has the following very compli¬ 
cated functional structure PS [35] [36l [37]: 


n MF 

t'ka 




oj — ( nt a E _ + n_ a E + ) — A (k) 

(u- E + - nI CT Ai(fc))(w - E_- n+ CT A 2 (&)) - nZ (T nt a \‘z{k)\±{k) 


(78) 


Here, the quantities A i(k) are the components of the generalized mean field, which cannot be 
reduced to the functional of the mean particle’s densities. The expression for Green function (1781) 
can be written down in the form of the following generalized two-pole solution 


n MF 


M 


nt a (l + cb x ) n_ a (l + da x ) 
a — db~ 1 c b — ca~ 1 d 


n_ 


u-E_- nZ a W k _ a 


+ 


n_ 


u- E+- n_„W, 


—o' vv k—o 


(79) 
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where 


nt (T n_ (J W±_ a = N 1 ^Uj exp [~ik{Ri - Rj)\ x (80) 

U 

(((4_ CT n^Oj_ CT ) + (a^n? a]_ ff )) + 

Green function (1791) is the most general solution of the Hubbard model within the generalized mean 
field approximation. Equation (1801) is nothing else but the explicit expression for the generalized 
mean field. As we see, this mean field is not a functional of the mean particle’s densities. The 
solution (1791) is more general than the solution ” Hubbard IIP m and other two-pole solutions. 
Hence it was shown in the papers [171 35l 1361137| that the solution ” Hubbard i” m is a particular 
case of the solution (1791) . which corresponds to the additional approximation 


n± a n_ a W ± (k ) ~ N 1 ^ tijesxp[-ik(Ri - Rj)]{nf_ a nf_ ff ). (81) 

ij 


Assuming ( nj- a ni- a ) ~ n 2 _ a , we obtain the approximation ” Hubbard p da Ea eh ez!. Thus, we 
have shown that in the cases of systems of strongly correlated particles with a complicated char¬ 
acter of quasiparticle spectrums the generalized mean fields can have quite a nontrivial structure, 
which is difficult to establish by using any kind of independent considerations. The method of 
irreducible Green functions allows one to obtain this structure in the most general form. 

One should note that the BCS-Bogoliubov superconductivity theory pia lama da is formulated 
in terms of a trial (approximating) Hamiltonian H mo d, which is a quadratic form with respect 
to the second-quantized creation and annihilation operators, including the terms responsible for 
anomalous (or non-diagonal) averages. For the single-band Hubbard model the BCS-Bogoliubov 
functional of generalized mean fields can be written in the following form [92 ,, 93; [Ml [95] 


T, C =U 


(®j — a Qi—cr) {dia&i—cr) 

-<«L4r> -<°L a i«) , 


(82) 


The anomalous (or nondiagonal) mean values in this expression fix the vacuum state of the system 
exactly in the BCS-Bogoliubov form. 

It is worth mentioning that the modern microscopic theory of superconductivity was given a 
rigorous mathematical formulation in the classic works of Bogoliubov and co-workers PEI El PH 
CEO simultaneously with the Bardeen, Cooper, and Sehrieffer (BCS) theory. It was shown that the 
equations of superconductivity can be derived from the fundamental electron-ion and electron- 
electron interactions. The set of equations obtained is known as the Eliashberg equations. They 
enable us to investigate the electronic and lattice properties of a metal in both the normal and 
superconducting states. Moreover, the Eliashberg equations are appropriate to the description 
of strong coupling superconductors, in contrast to the equations, which are valid in the weak 
coupling regime and describe the electron subsystem in the superconducting state only. 

In paper [92J on the basis of the BCS-Bogoliubov functional of generalized mean fields a system 
of equations of superconductivity for the tight-binding electrons in the transition metal described 
by the Hubbard Hamiltonian was derived. The electron-phonon interaction was written down for 
the ’’rigid ion model”. Neglecting the vertex corrections in the self-energy operator the closed 
system of equations was obtained. 

In paper [93| this approach was extended for the Barisic-Labbe-Friedel model of a transition metal. 
The renormalized electron and phonon spectra of the model were derived using the method of 
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irreducible Green functions pH m EH ED in a the self-consistent way. For the band and atomic 
limits of the Hubbard model the explicit solutions for the electron and phonon energies were 
obtained. The energy gap, appearing between electron bands in the strong correlation limit, 
persists in that calculations. The Eliashberg-type equations of superconductivity were obtained 
also. 

The equations of strong coupling superconductivity in disordered transition metal alloys have 
been derived in paper [94] by means of irreducible Green functions method and on the basis of the 
alloy version of the Barisic-Labbe-Friedel model for electron-ion interaction. The configurational 
averaging has been performed by means of the coherent potential approximation. Making some 
approximations, the formulas for the superconducting transition temperature Tq and the electron- 
phonon coupling constant have been obtained. These depend on the alloy component and total 
densities of states, the phonon Green function, and the parameters of the model. 

To summarize, various schemes of ’’effective mean-field theory” taking into account correlations 
were proposed [35j [36l EH [96l EZUMJ El [Ml Ml Ml Ml Ml Ml M] • The main efforts were 
directed to the aim to describe suitably the collective behavior of particles in terms of effective- 
field distribution which satisfies a self-consistent condition. However, although the self-consistent 
field approximation often is a reasonable approximation away from the critical point, it usually 
breaks down in its immediate neighborhood. 

It is of importance to stress again that from our point of view, in real mean-field theory, the mean 
field appearing in the single-site problem should be a scalar or vectorial time-independent 
quantity. 

6 Symmetry Broken Solutions 

The formalism of the previous sections may be extended to incorporate the broken symmetry 
solutions pT, 18j 107 of the interacting many particle systems, e.g. the pairing effects present in 
superconductors [a hi ns], etc. Our purpose in this section is to attract attention to subtle points 
which are essential for establishing a connection of the generalized mean field approximation and 
the broken symmetry solutions [171 fl8 . 10 ? . 

It is well known that a symmetry can be exact or approximate. Symmetries inherent in the phys¬ 
ical laws may be dynamically and spontaneously broken, i.e., they may not manifest themselves 
in the actual phenomena. It can be as well broken by certain reasons [108 , 109 ]. 

Within the literature the term broken symmetry is used both very often and with different mean¬ 
ings. There are two terms, the spontaneous breakdown of symmetries and dynamical symmetry 
breaking, which sometimes have been used as opposed but such a distinction is irrelevant. How¬ 
ever, the two terms may be used interchangeably. It should be stressed that a symmetry implies 
degeneracy. In general there are a multiplets of equivalent states related to each other by congru¬ 
ence operations. They can be distinguished only relative to a weakly coupled external environment 
which breaks the symmetry. Local gauged symmetries, however, cannot be broken this way be¬ 
cause such an extended environment is not allowed (a superselection rule), so all states are singlets, 
i.e., the multiplicities are not observable except possibly for their global part. 

It is known that when the Hamiltonian of a system is invariant under a symmetry operation, but 
the ground state is not, the symmetry of the system can be spontaneously broken. Symmetry 
breaking is termed spontaneous when there is no explicit term in a Lagrangian which manifestly 
breaks the symmetry. 

Peierls mm gave a general definition of the notion of the spontaneous breakdown of sym¬ 
metries which is suited equally well for the physics of particles and condensed matter physics. 
According to Peierls mol eh, the term broken symmetries relates to situations in which symme¬ 
tries which we expect to hold are valid only approximately or fail completely in certain situations. 
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The intriguing mechanism of spontaneous symmetry breaking is a unifying concept that lie at 
the basis of most of the recent developments in theoretical physics, from statistical mechanics to 
many-body theory and to elementary particles theory [1081 [109) . The existence of degeneracy 
in the energy states of a quantal system is related to the invariance or symmetry properties of 
the system. By applying the symmetry operation to the ground state, one can transform it to a 
different but equivalent ground state. Thus the ground state is degenerate, and in the case of a 
continuous symmetry, infinitely degenerate. The real, or relevant, ground state of the system can 
only be one of these degenerate states. A system may exhibit the full symmetry of its Lagrangian, 
but it is characteristic of infinitely large systems that they also may condense into states of lower 
symmetry. 

It should be pointed out that Bogoliubov’s method of quasiaverages mm [15] gives the deep 
foundation and clarification of the concept of broken symmetry. It makes the emphasis on the 
notion of degeneracy and plays an important role in equilibrium statistical mechanics of many- 
particle systems. According to that concept, infinitely small perturbations can trigger macroscopic 
responses in the system if they break some symmetry and remove the related degeneracy (or quasi¬ 
degeneracy) of the equilibrium state. As a result, they can produce macroscopic effects even when 
the perturbation magnitude is tend to zero, provided that happens after passing to the ther¬ 
modynamic limit [45] . This approach has penetrated, directly or indirectly, many areas of the 
contemporary physics. 

The article [18] examines the Bogoliubov’s notion of quasiaverages, from the original papers [1], 
through to modern theoretical concepts and ideas of how to describe both the degeneracy, broken 
symmetry and the diversity of the energy scales in the many-particle interacting systems. Current 
trends for extending and using Bogoliubov’s ideas to quantum field theory and condensed matter 
physics problems were discussed, including microscopic theory of superfluidity and superconduc¬ 
tivity, quantum theory of magnetism of complex materials, Bose-Einstein condensation, chirality 
of molecules, etc. Practical techniques covered include quasiaverages, Bogoliubov theorem on the 
singularity of 1 /q 2 , Bogoliubov’s inequality, and its applications to condensed matter physics. 

It was demonstrated there that the profound and innovative idea of quasiaverages formulated by 
N.N. Bogoliubov, gives the so-called macro-objectivation of the degeneracy in domain of quantum 
statistical mechanics, quantum field theory and in the quantum physics in general. 

The quasiaverages may be obtained from the ordinary averages by using the cluster property 
which was formulated by Bogoliubov MES]. This was first done when deriving the Boltz¬ 
mann equations from the chain of equations for distribution functions, and in the investigation of 
the model Hamiltonian in the theory of superconductivity [SUES]- To demonstrate this let us 
consider averages (quasiaverages) of the form 

F(ti,xi,...t n ,x n ) = (... #(ii,xi)... ^(tj,Xj )...), (83) 

where the number of creation operators may be not equal to the number of annihilation opera¬ 
tors T. We fix times and split the arguments (t\, x\,... t n , x n ) into several clusters (..., t a ,x a ,...),..., 
(..., tp, xp ,...). Then it is reasonably to assume that the distances between all clusters \x a — xp\ 
tend to infinity. Then, according to the cluster property, the average value (1831) tends to the prod¬ 
uct of averages of collections of operators with the arguments (..., t a ,x a ,...),..., (..., tp, xp ,...) 

lirn F(ti,xi,... t n ,x n ) = F(. ..,t a ,x a ,...)... F(...,tp,xp ,...). (84) 

\X a —Xp\-¥<X 

For equilibrium states with small densities and short-range potential, the validity of this property 
can be proved mmm- For the general case, the validity of the cluster property has not yet 
been proved. Bogoliubov formulated it not only for ordinary averages but also for quasiaverages, 
i.e., for anomalous averages, too. It works for many important models, including the models of 
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superfluidity and superconductivity la m la¬ 
in his work The Theory of Superfluidity m, Bogoliubov gave a microscopic explanation of the 
phenomenon of superfluidity [2l 1113] . Before his works, there were the phenomenological theories 
which were based on an assumption about the form of the spectrum of elementary excitations. 
Bogoliubov has started from the general Hamiltonian for Bose systems and assumed that a macro¬ 
scopic number of particles are found in the ground state with zero momentum, and therefore the 
creation and annihilation operators of particles with zero momentum are c-numbers m- As a 
result a definite approximating Hamiltonian was obtained, consisting from a quadratic form of the 
creation and annihilation operators. The usual perturbation theory proved to be inapplicable to it 
because of the strong interaction of particles with opposite momenta. Therefore the Hamiltonian 
was diagonalized with the help of the canonical transformations (the Bogolyubov u — v transfor¬ 
mations). This permitted one to calculate the spectrum of elementary perturbations outside the 
framework of perturbation theory. Decomposing the field operators into c-numerical and operator 
parts, Bogoliubov in fact introduced into quantum theory the method of spontaneous symmetry 
breakdown for systems with degenerate ground state. This method was rediscovered in quantum 
field theory a decade later [T8| . 

To illustrate these statements consider Bogoliubov’s theory of a Bose-system with separated con¬ 
densate, which is given by the Hamiltonian 13 SI US 

H\= [ T^(x)(-)T(x)dx —/r [ (85) 

Ja 2 m J A 

+7: f ^ t (xi)^ t (x 2 )$(xi - x 2 )Af{x 2 )^(x 1 )dx 1 dx 2 . 

2 J\2 

This Hamiltonian can be written also in the following form 

H A = H 0 +Hi = j^(q)(—^-)^>(q)dq ( 86 ) 

+77 [ ^\q)^{q')^(q-q')^W)^(q)dqdq'. 

2 J A2 

Here, it (<7), and !l( q ) are the operators of annihilation and creation of bosons. They satisfy the 
canonical commutation relations 

[*(<?), ^W)} = % - </); [*(?), *(</)] = [* f (?), *V)] = 0 . (87) 

The system of bosons is contained in the cube A with the edge L and volume V . It was assumed 
that it satisfies periodic boundary conditions and the potential &(q) is spherically symmetric and 
proportional to the small parameter. It was also assumed that, at temperature zero, a certain 
macroscopic number of particles having a nonzero density is situated in the state with momentum 
zero. 

The operators T( q ), and ^'(q) are represented in the form 

*(q) = “o/W; tf(q) = 4/W, (88) 

where ao and aj are the operators of annihilation and creation of particles with momentum zero. 
To explain the phenomenon of superfluidity [3 1112] . one should calculate the spectrum of the 
Hamiltonian, which is quite a difficult problem. Bogoliubov suggested the idea of approximate 
calculation of the spectrum of the ground state and its elementary excitations based on the physical 
nature of superfluidity. His idea consists of a few assumptions. The main assumption is that at 
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temperature zero the macroscopic number of particles (with nonzero density) has the momentum 
zero. Therefore, in the thermodynamic limit [35J, the operators a^/y/V and Oq/'JV commute 


linr 

V—>-oo 


a 0 /W, al/y/V 



(89) 


and are c-nunrbers. Hence, the operator of the number of particles N ,o = a^ao is a c-nunrber, too. 
D. Ya. Petrina [115] shed an additional light on the problem of an approximation of general 
Hamiltonians by Hamiltonians of the theories of superconductivity and superfluidity. In his highly 
interesting paper [115] Petrina pointed out that the model Hamiltonian of the theory of supercon¬ 
ductivity mm can be obtained from the general Hamiltonian for Fermi systems if the Kronecker 
symbol, which expresses the law of conservation of momentum in the interaction Hamiltonian, is 
replaced by two Kronecker symbols so that only particles with opposite momenta interact. The 
model Hamiltonian of the theory of superfluidity can be obtained from the general Hamiltonian 
for Bose systems if we replace the Kronecker symbol, which expresses the law of conservation of 
momentum, by several Kronecker symbols, preserving only the terms that contain at least two 
operators with momenta zero in the interaction Hamiltonian. This list of model systems can be 
continued m- 

The concept of quasiaverages was introduced by Bogoliubov on the basis of an analysis of rnany- 
particle systems with a degenerate statistical equilibrium state. Such states are inherent to various 
physical many-particle systems. Those are liquid helium in the superfluid phase, metals in the 
superconducting state, magnets in the ferromagnetically ordered state, liquid crystal states, the 
states of superfluid nuclear matter, etc. 

In many-body interacting systems, the symmetry is important in classifying different phases and 
in understanding the phase transitions between them. According to Bogoliubov’s ideas mm 
B> 112, 1107 1 in each condensed phase, in addition to the normal process, there is an anomalous 
process (or processes) which can take place because of the long-range internal field, with a corre¬ 
sponding propagator. Additionally, the Goldstone theorem [18] states that, in a system in which 
a continuous symmetry is broken ( i.e. a system such that the ground state is not invariant under 
the operations of a continuous unitary group whose generators commute with the Hamiltonian), 
there exists a collective mode with frequency vanishing, as the momentum goes to zero. For many- 
particle systems on a lattice, this statement needs a proper adaptation. In the above form, the 
Goldstone theorem is true only if the condensed and normal phases have the same translational 
properties. When translational symmetry is also broken, the Goldstone mode appears at a zero 
frequency but at nonzero momentum, e.g., a crystal and a helical spin-density-wave ordering (see 
for discussion Refs. puma). 

The antiferromagnetic state is characterized by a spatially changing component of magnetization 
which varies in such a way that the net magnetization of the system is zero. The concept of anti- 
ferromagnetism of localized spins which is based on the Heisenberg model and the two-sublattice 
Neel ground state is relatively well founded contrary to the antiferromagnetism of delocalized 
or itinerant electrons. The itinerant-electron picture is the alternative conceptual picture for 
magnetism [118] , In the antiferromagnetic many-body problem there is an additional ’’symme¬ 
try broken” aspect mm- For an antiferromagnet, contrary to ferromagnet, the one-electron 
Hartree-Fock potential can violate the translational crystal symmetry. The period of the anti¬ 
ferromagnetic spin structure L is greater than the lattice constant a. The Hartree-Fock is the 
simplest approximation but neglects the important dynamical part. To include the dynamics one 
should take into consideration the correlation effects. 

The anomalous propagators for an interacting many-fermion system corresponding to the ferro¬ 
magnetic (FM) , antiferromagnetic (AFM), and superconducting (SC) long-range ordering are 
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given by 


FM : Gf m ~ ((a ka ; a\_ a )), (90) 

AFM : G a f m ~ ((fl-fc+Qo-j ^k+Q'a'^ ’ 

SC : G sc ~ (( aka\a-k-a))■ 


In the spin-density-wave case, a particle picks up a momentum Q — Q' from scattering against 
the periodic structure of the spiral (nonuniform) internal field, and has its spin changed from a to 
cr' by the spin-aligning character of the internal field. The Long-Range-Order (LRO) parameters 
are: 


FM : m = l/N^(a\ a a k - a ), (91) 

kcr 

AFM : Mq = ^2(al a a k+Q _ a ), 

kcr 

SC: A = J>!. H 4 t }. 

k 

It is of importance to note that the long-range order parameters are functions of the internal field, 
which is itself a function of the order parameter. There is a more mathematical way of formulating 
this assertion. As it was stressed earlier EH), the notion symmetry breaking means that the state 
fails to have the symmetry that the Hamiltonian has. 

In terms of the theory of quasiaverages, a true breaking of symmetry can arise only if there are in¬ 
finitesimal ’’source fields”. Indeed, for the rotationally and translationally invariant Hamiltonian, 
suitable source terms should be added: 

FM : £HbH x y j al^ak-g, (92) 

kcr 

AFM : £^ B H^a\ a a k+Q _ a , 

kQ 

SC : ev E 

k 

where £ —>• 0 is to be taken at the end of calculations. 

For example, broken symmetry solutions of the spin-density-wave type imply that the vector Q is 
a measure of the inhomogeneity or breaking of translational symmetry. 

In this context the Hubbard model is a very interesting tool for analyzing the broken symmetry 
concept [35lI361I3T] . It is possible to show that antiferromagnetic state and more complicated states 
(e.g. ferrimagnetic) can be made eigenfunctions of the self-consistent field equations within an 
’’extended” (or generalized) mean-field approach, assuming that the anomalous averages (a| CT aj_ CT ) 
determine the behaviour of the system on the same footing as the ’’normal” density of quasipar¬ 
ticles (aj Oj CT ). It is clear, however, that these “spin-flip” terms break the rotational symmetry of 
the Hubbard Hamiltonian. For the single-band Hubbard Hamiltonian, the averages (aj aj i(T ) = 0 
because of the rotational symmetry of the Hubbard model. The inclusion of anomalous averages 
leads to the following approximation 

^i—cr^icy ~ (^z—cr)^icr (Q i i_ (7 Q j ia')Q i i—cr‘ (^^) 
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Thus, in addition to the standard Hartree-Fock term, the new so-called “spin-flip” terms are 
retained m- This example clearly shows that the structure of mean field follows from the 
specificity of the problem and should be defined in a proper way. So, one needs a properly 
defined effective Hamiltonian H e R. In paper [117] we thoroughly analyzed the proper definition 
of the irreducible Green functions which includes the “spin-flip” terms for the case of itinerant 
antiferromagnetism of correlated lattice fermions. For the single-orbital Hubbard model [35l L36j 
mm, the definition of the ’’irreducible” part should be modified in the following way: 

\{ a k+po a \,+q-a a q-cr\ a \ a ))u = {{ a k+pa a \>+q- a a q-<y\ a \o))u ~~ 

fip,o{W'q—cr')C'ka {^k+pa^p-\-q—a }(( a q-cr| a [. cr })w (94) 

From this definition it follows that this way of introduction of the irreducible Green functions 
broadens the initial algebra of operators and the initial set of the Green functions. This means 
that the “actual” algebra of operators must include the spin-flip terms from the beginning, namely: 
{ai a , aj CT , nia, aj^aj-o-). The corresponding initial Green function will be of the form 

V,>, (0 ' 0®) 

With this definition, one introduces the so-called anomalous (off-diagonal) Green functions which 
fix the relevant vacuum and select the proper symmetry broken solutions. In fact, this approxi¬ 
mation was investigated earlier by Kishore and Joshi [ 119] . They clearly pointed out that they 
assumed a system to be magnetized in the x direction instead of the conventional 2 axis. 

The problem of finding the superconducting, ferromagnetic and antiferromagnetic ’’symmetry bro¬ 
ken” solutions of the correlated lattice fermion models within irreducible Green functions method 
was investigated in Refs. [17, [35J S5, 3711117] . A unified scheme for the construction of gener¬ 
alized mean fields (elastic scattering corrections) and self-energy (inelastic scattering) in terms 
of the Dyson equation was generalized in order to include the ’’source fields”. The ’’symmetry 
broken” dynamic solutions of the Hubbard model which correspond to various types of itinerant 
antiferromagnetism were discussed as well [lTl [H [36l [37,, 117] . This approach complements previ¬ 
ous studies of microscopic theory of the Heisenberg antiferromagnet [120] and clarifies the concepts 
of Neel sublattices for localized and itinerant antiferromagnetism and ’’spin-aligning fields” of cor¬ 
related lattice fermions. 

We shall see shortly that in order to discuss the mean field theory (and generalized mean fields) on 
the firm ground the Bogoliubov inequality provides the formal basis and effective general approach. 

7 The Mathematical Tools 

Before entering fully into our subject, we must recall some basic statements. This will be necessary 
for the following discussion. 

The number of inequalities in mathematical physics is extraordinary plentiful and the literature on 
inequalities is vast p2H EH [Ml EH ESS EH EH EH QH EH EH EH m • The physicists 
are interested mostly in intuitive, physical forms of inequalities rather than in their most general 
versions. Often it is easier to catch the beauty and importance of original versions rather than 
decoding their later, abstract forms. 

Many inequalities are of a great use and directly related with the notion of entropy, especially 
with quantum entropy [12411134] , The von Neuman entropy of p € S n , S(p), is defined by 

S(p) = -Tr(plogp). (96) 
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The operator p log p is defined using the spectral theorem |124| . Here S n denotes the set of density 
matrices p on C n . 

In fact, S(p) depends on p only through its eigenvalues. 

n 

S(p) = -J2 XjiogXj. (97) 

3 = 1 

Otherwise put, the von Neumann entropy is unitarily invariant; i.e., 

S{UpU*) = S(p). (98) 

The convexity condition leads to m 

- S(p) = -log(ra). (99) 

This equality valid if and only if each A j = 1/n Thus, one may arrive at [ 124J 

0 < S{p) < logn ( 100 ) 


for all p 6 S n , and there is equality on the left if p is a pure state, and there is equality on the 
right if p = ( 1 /n)/. Actually, S(p) is not only a strictly concave function of the eigenvalues of p , 
it is strictly concave function of p itself. 

The notions of convexity and concavity of trace functions [124] are of great importance in math¬ 
ematical physics [S3, Gill- Inequalities for quantum mechanical entropies and related concave 
trace functions play a fundamental role in quantum information theory as well [124 . 134 j. 

A function / is convex in a given interval if its second derivative is always of the same sign in that 
interval. The sign of the second derivative can be chosen as positive (by multiplying by (—1) if 
necessary). Indeed, the notion of convexity means that if d 2 f/dx 2 > 0 in a given interval, Xj are 
a set of points in that interval, pj are a set of weights such that pj > 0 , which have the property 
T,jPj = 1 , then 

^Pjf{xj) > f(j2 p i x i)- ( 101 ) 

3 3 

The equality will be valid only if Xj = (x) = YljPj x j■ I n other words, a real-valued function 
f(x) defined on an interval is called convex (or convex downward or concave upward) if the line 
segment between any two points on the graph of the function lies above the graph, in a Euclidean 
space (or more generally a vector space) of at least two dimensions. Equivalently, a function is 
convex if its epigraph (the set of points on or above the graph of the function) is a convex set. 

A real-valued function / on an interval (or, more generally, a convex set in vector space) is said 
to be concave if, for any x\ and x 2 in the interval and for any a in [0,1], 

/((1 - a)xi + (a)x 2 ) > (1 - a)f(x i) + (a)f(x 2 ). (102) 


A function f(x) is concave over a convex set if and only if the function —f(x) is a convex function 
over the set. 

As an example we mentioned above briefly a reason this concavity matters, pointing to the in¬ 
equality (HOOD that was deduced from the concavity of the entropy as a function of the eigenvalues 
of p. 

It is of importance to stress that in quantum statistical mechanics, equilibrium states are deter¬ 
mined by maximum entropy principles [124] . and the fact that 


sup S (p) 


p£ S n 


log n, 


(103) 
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reflects the famous Boltzmann formulae 


S = k B \ogW. (104) 

It follows from Boltzmann definition that the entropy is larger if p is smeared out, where p is the 
probability density on phase space. The microscopic definition of entropy given by Boltzmann 
does not, by itself, explain the second law of thermodynamics, even in classical physics. The task 
to formulate these questions in a quantum framework was addressed by Oskar Klein in his seminal 
paper [ 137] of 1931. He found a fundamentally new way for information to be lost hence entropy 
to increase, special to quantum mechanics. This result was called Klein lemma (136111381 Il37j . 

M. B. Ruskai [138) has reviewed many fundamental properties of the quantum entropy [ 134] in¬ 
cluding one important class of inequalities relates the entropy of subsystems to that of a composite 
system. That article presented self-contained proofs of the strong subadditivity inequality for von 
Neumann quantum entropy, S(p), and some related inequalities for the quantum relative entropy, 
most notably its convexity and its monotonicity under stochastic maps. The approach to subad¬ 
ditivity and relative entropy presented was used to obtain conditions for equality in properties of 
relative entropy, including its joint convexity and monotonicity. In addition, the Klein inequality 
was presented there in detail. 

Indeed, the fact that the relative entropy is positive [ 138] . i.e., H(p\,p 2 ) > 0 when Tr/q = Tr/q, 
is an immediate consequence of the following fundamental convexity result due to Klein ]137 ., 139 . 
1140) . The corresponding theorem [ 138] states that for A, B > 0 

TrH(logH — log-B) > Tr(H — B), (105) 

with equality if and only if (A = B). 

In more general form )12'4j the Klein inequality may be formulated in the following way. 

For all A,B£ H„ , and all differentiable convex functions / : M — > K, or for all A, B £ and all 
differentiable convex functions / : (0, oo) —>• M 

Tr (/(H) - f(B) -(A- B)f'(B )) > 0. (106) 

In either case, if / is strictly convex, there is equality if and only if A = B. 

A few more words about Oskar Klein and his inequality will not be out of place here. Oskar 
Klein (1894 - 1977) was famous Swedish theoretical physicist [] who worked on a wide variety 
of subjects [141 ]. For example, the Klein-Gordon equation was the first relativistic wave equa¬ 
tion. Oskar Klein was also a collaborator of Niels Bohr in Copenhagen. It is interesting to note 
that Oskar Klein defended his thesis and was awarded his doctoral degree in 1921 for his work 
in physical chemistry about strong electrolytes. In 1931 Oskar Klein [137) 1391 140;. 141] . using 
his experience in both quantum and statistical mechanics, succeeded in solving the problem of 
whether the quantum statistics on molecular level can explain how the entropy increases with 
time in accordance with the second law of thermodynamics. The problem in classical statistical 
mechanics had been already noticed by Gibbs earlier. Klein proof [1371 Il39l 140) . that used the 
statement that only the diagonal elements in the density matrix for the phase space of the particles 
are relevant for the entropy, has led him to the Klein lemma. With Klein lemma, the entropy can 
increase according to the formula of Boltzmann microscopic definition, where it is described with 
the number of states in the phase space. A useful and informative discussion of the Klein paper 
and Klein lemma was carried out in the book of R. Jancel |il42j . 

According to M. B. Ruskai [138] . the closely related Peierls-Bogoliubov inequality is sometimes 
used instead of Klein inequality. Golden-Thompson and Peierls-Bogoliubov inequalities were ex¬ 
tended to von Neumann algebras, which have traces, by Ruskai [ 128 ] (see also Ref. [ 143] ). H. 
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Araki [ 129] extended them to a general von Neumann algebra. This kind of investigations is 
particularly valuable since the Bogoliubov inequality is remarkable because it have significant 
applications in statistical quantum mechanics mmmmam- It provides insight into a 
number of other interesting questions as well. 

It will be of use to write down the mathematical formulation of Peierls-Bogoliubov inequalities 
which was provided by Carlen [ 124] , 

Let us consider A € H n , and let / be any convex function on M. Let {u\,... ,u n } be any or¬ 
thonormal base of C n . Then 


< Tr[/(A)]. (107) 

3 =1 

There is equality if each Uj is an eigenvector of A, and if / is strictly convex, only in this case. 
Now consider the formulation of the generalized Peierls-Bogoliubov inequality [124] . For every 
natural number n, the map A H» log^Tr[exp(A)]^ is convex on H n . As a consequence one may 
deduce [124] that 


/Tr[exp(A + B)] \ Tr[Hexp(A)] 
° v Tr[exp(A)] / — Tr[exp(A)] 


(108) 


Frequently this relation, which has many uses, is referred to as the Peierls-Bogoliubov inequality. 
It is worth noting that according to tradition the term Gibbs-Bogoliubov inequality [8] is used for a 
classical statistical mechanical systems and term Peierls-Bogoliubov inequality [124] for quantum 
statistical mechanical systems. 

At the very least, it must have been meant to indicate that Peierls inequality does not have a 
classical analog, whereas Bogoliubov inequality has. 


8 Variational Principle of Bogoliubov 

It is known that there are several variational principles which provide upper bounds for the 
Helmholtz free-energy function. With these instruments, it is possible to construct various approx¬ 
imations to the statistical thermodynamic behavior of systems. For any variational formulation, 
its effectiveness as a minimal principle will be enhanced considerably if there is a workable tool 
for determining lower bounds to the Helmholtz free energy function. Bogoliubov inequality for 
the free-energy functional is an inequality that gives rise to a variational principle of statistical 
mechanics. It is used [H13 El SI E] to obtain the exact thermodynamic limit [IS] solutions of 
model problems in statistical physics, in studies using the method of molecular fields, in proving 
the existence of the thermodynamic limit [35], and also in order to obtain physically important 
estimates for the free energies of various many-particle interacting systems. 

A clear formulation of the variational principle of Bogoliubov and Bogoliubov inequality for the 
free-energy functional was carried out by S. V. Tyablikov [5j. We shall follow close to that formu¬ 
lation. Tyablikov [5j used the theorems relating to the minimum values of the free energy. As a 
result, it was possible to formulate a principle which then was used to deduce the molecular field 
equations. 

Principle of the free energy minimum is based on the following arguments. Let us consider an 
arbitrary complete system of orthonormalized functions which are not the eigenfunctions 

of the Hamiltonian Ti of a system. Then it is possible to write down the inequality 

F(H) < F mod (n). (109) 

Here F(H) is the intrinsic free energy of the system 

F(M) = -einZ, Z = J^exp (~E u /d), (110) 
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0 = ksT, E u are eigenfunctions of the Hamiltonian 77, F mod {T-L) is the model free energy, which 
gives approximately the upper limit of the intrinsic free energy: 


T m od(77) = -01nZ mod , Z mod = ^2exp(-'H nn /6), U nn = (^*,77</j n ). 

n 

The inequality (I109D may be written also in the following way: 

Z > ^mod- 


( 111 ) 


( 112 ) 


The relationships represented by the equality sign in Eqs. (11091) and (I112D applies if ip n are eigen¬ 
functions of the Hamiltonian of the system. It should be noted that for finite values of the number 
of partial sums Z^ N \ the quantity F^ d does not reach its maximum for any system of functions 
ipi,... ,ipN- In fact, the inequality will be satisfied really BM in the limit N —» oo. 

Using these results, it is possible to formulate a variational principle for the approximate deter¬ 
mination of the free energy of a system [5j. To proceed, let us suppose that the functions {tyn) 
depend on some arbitrary parameter A. It was established above that 

F(H) < F mod (H) = -01n^exp(-77 nri (A)/0). (113) 

n 

It is clear that the best approximation for the upper limit of the free energy F is obtained by 
selecting the values of the parameter A in accordance with the condition for the minimum of the 
model free energy F mod . Indeed, let the Hamiltonian of the system, 77, be written in the form 

77 = 77 0 (A) + A77(A) = 77 0 (A) + (n- 74(A)) , (114) 

where 74(A) is some operator depending on the parameter A. The concrete form of the operator 
74(A) should be selected on the basis of convenience in calculations. We shall use notation E® 
and ip n for the eigenvalues and the eigenfunctions of the operator 74- To denote the diagonal 
matrix elements of the operator A4 in terms of the functions ip n we shall use the notation A74n- 
For a generality, we shall assume that <p n are not the eigenfunctions of the total Hamiltonian 77. 
Clearly, E® and A77 nn are also some functions of the parameter A. In this sense, the system of 
functions {^Pn} plays a role of a trial system of functions. Then we may write that 

H nn = E° n + A74n = E° n + (74n - £°) • (115) 

As a consequence, the free energy will satisfy the inequality 

F{H) < -01nJ]exp-(E° + A74„)i (116) 


Now let us suppose that the operator AT l can be considered as a small perturbation compared 
with the operator 77. We obtain then [5j, to within quantities of the first order of smallness with 
respect to A77 


F(H) < F(77 0 ) + 


Tr( A77 exp(—77 q 


(117) 


Trl exp(—77 q 


Note that in this case, the best approximation to the upper limit of the free energy is obtained 
by selecting the value of the parameter A from the condition for the minimum of the right-hand 
side of Eq. (I117D . The formulation of the variational principle of Eq. (I117D is more restricted than 
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the initial formulation of Eq. (110911 . 

The variational principle in the form of Eq. (1117f) can be strengthened, following to the Bogoliubov 
suggestion |5], by removing the limitation of the smallness of the operator AH. As a result we 
obtain 

F(H) < F mod (H). 


Here 


^mod(W) = F(H 0 ) + 


Tr (AH exp(— Hq 


Tr^exp (~H 0 /6) 

F(Hq) = —6lnTrexp(—'Ho/9). 

Hence one may write down also that for a system with the Hamiltonian 

n = n 0 + ah 

the free energy has a certain upper bound. Bogoliubov inequality states that: 


(118) 

(119) 

( 120 ) 

( 121 ) 


F<F 0 + (H-Ho)o, (122) 

or 

F<F q (H)o-TSo, (123) 

where Sq is the entropy and where the average is taken over the equilibrium ensemble of the 
reference system with Hamiltonian Hq. Usually "Ho contains one or more variational parameters 
which are chosen such as to minimize the right hand side of Eci. (H22D . In the special case that 
the reference Hamiltonian is that of a non-interacting system and can thus be written as a sum 
single-particles Hamiltonians [5] 

N 

-Ho = J2 h i- ( 124 ) 

i —1 

Then it is possible to improve the upper bound by minimizing the right hand side of the inequality 
(I122D . The minimizing reference system is then the trial approximation to the true system using 
non-correlated degrees of freedom, and is known as the mean field approximation. 

Starting with the one-particle model Hamiltonian that can be exactly solved in the Bogoliubov 
variational method, one may get a self-consistent result such as the molecular field theory in the 
ferromagnet and the Hartree-Fock approximation in many-particle problems. Since the variational 
method yields a result which is always greater than the correct answer, the mathematical meaning 
for improving upon the approximation in the variational method is strictly defined by lowering 
the upper bound of the free energy. But these variational methods, the molecular field theory 
and the Hartree-Fock approximation, have such a feature that the correlation effects cannot be 
taken into account correctly. In general case [5] the Hamiltonian of a system contains interpar¬ 
ticle interactions. Thus Bogoliubov variational principle can be considered as the mathematical 
foundation of the mean field approximation in the theory of many-particle interacting systems. 
Using the Klein inequality (110611 it is possible to write down a general form of the Bogoliubov 
inequality for the free energy functional. The following inequality is valid for any Hermitian 
operators and Hi and H 2 

N-^Hi - H 2 ) Hl < (/(Hi) - /(H 2 )) < N-\Hi - H 2 ) H2 , (125) 

where 

f(H) = —6W _1 lnTrexp(—H/0). (126) 
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This expression has the meaning of the free energy density for a system with Hamiltonian H and 
the extensive parameter N may be treated as the number of particles or the volume, depending 
on the system. 

Derrick [ 146] established a simple variational bound to the entropy S(E ) of a system with energy 


E 


S(E) > -Jfe B ln(Trl/ 2 ) 


(127) 


for all Hermitian matrices U (with no negative eigenvalues) for which Tr?7 = 1 and Tr(HU) = E, 
where H is the Hamiltonian. 

This principle has the advantage that U 2 is in general much easier to evaluate than U In U which 
appears in the conventional bound given by von Neumann 


S(E) > -k B (TvU In £/). (128) 

There are numerous methods for proving of the Bogoliubov inequality [5;, [Til, HI LE2J [Ml IE1 d47 „ 
1148] , A. Oguchi [149] proposed an approach for determination of an upper bound and a lower 
bound of the Helmholtz free energy in the statistical physics. He used as a basic tool the Klein 
lemma [1244138 . 137 ]. He obtained a new approximate expression of the free energy. This approx¬ 
imate value of the free energy was conjectured to be greater than the lower bound and less than 
the upper bound. An approach which can be extended to improve the approximation was formu¬ 
lated. The upper bound and the lower bound of the approximate free energy converge to the true 
free energy as the successive approximation proceeds. The method was first applied to the Ising 
ferromagnet and then applied to the Heisenberg ferromagnet. In the simplest approximation the 
results agrees with the Bethe-Peierls approximation for the Ising model and the constant coupling 
approximation for the Heisenberg model. In his subsequent paper, A. Oguchi [150] formulated a 
new variational method for the free energy in statistical physics. According to his calculations, the 
value of the free energy was obtained by using this new variational method was lower than that 
of the Bogolyubov variational method. Author concluded that the new variational free energy 
satisfies the thermodynamic stability criterion. 

However, J. Stolze [ 151] by careful examination of the papers [1491 1150] , has shown that the cal¬ 
culation in Ref. [ 150] contains a mistake which invalidates the result. He also pointed out several 
errors seriously affecting the results of an earlier paper [149] . Oguchi assumed that the Hamilto¬ 
nian T~Lo contains a variational parameter ”a” distributed according to a probability density P(a). 
Stolze derived a corrected inequality which clearly states that the new upper bound on the free 
energy suggested by Oguchi }l49l 1501] cannot be better (i.e., lower) than the Peierls-Bogoliubov 
bound, no matter how cleverly P(a ) was chosen. This shows clearly that no advantage over the 
Peierls-Bogolyubov bound was obtained. 

The standard proof was given in Callen second edition book on thermodynamics [48] for the 
case when the unperturbed Hamiltonian and the perturbation commute. Another proof (for the 
general case), was carried out in Feynman book on statistical mechanics [147] . Feynman used 
Baker-Campbell-Hausdorff expansion [123[ 1124] for the exponential of a sum of two non commut¬ 
ing operators. Prato and Barraco [148] presented a proof of the Bogoliubov inequality that does 
not require of the Baker-Campbell-Hausdorff expansion. 

Several variational approaches for the free energy have been proposed [ 152(1153 ] as attempts to im¬ 
prove results obtained through the well established Bogoliubov principle. This principle requires 
the use of a trial Hamiltonian depending on one or more variational parameters. The only way to 
improve the Bogoliubov principle by itself is to choose a more complete trial Hamiltonian, closing 
it to the exact one, but in almost all cases the possibilities are soon exhausted. The usual mean 
field approximation may be obtained using the above principle utilizing a sum of single spins in 
an effective held (the variational parameter) as the trial Hamiltonian. 
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Lowdin [ 154] and Lowdin and Nagel [155] studied a generalization of the Gibbs-Bogoliubov in¬ 
equality F < Fo + (T~L — 7io)o for the free energy F which leads to a variation principle for this 
quantity that may be of importance in certain computational applications to quantum systems. 
This approach is coupled with a study of the perturbation expansion of the free energy for a 
canonical ensemble with = T~Lq + XV in the general case when T~Lq and V do not commute. A 
simple proof was given for the thermodynamic inequality F — Fq — (7i — JLo)o < 0 in the case when 
the two Hamiltonian "Ho and V do not commute. The second- and high-order derivatives of the 
free energy with respect to the perturbation parameter A were calculated. From the second-order 
term was finally obtained a second-order correction to the previous variational minimum for the 
free energy. 

A. Decoster [39] established a sequence of inequalities which generalize the the Gibbs-Bogoliubov 
inequality in classical statistical mechanics and the Peierls and Bogoliubov inequalities in quantum 
mechanics; they can be presented as rearrangements of perturbation expansions, which provide 
exact bounds which are used in variational calculations. 

W. Kramarczyk [1561 argued that the Bogoliubov variational principle may be shown to be equiva¬ 
lent to the minimizing of the information gained while replacing the exact state by an approximate 
one. Consequently, the quasiparticles introduced in the thermal Hartree-Fock approximation may 
be redefined information-theoretically. 

9 Applications of the Bogoliubov Variational Principle 

Bogoliubov variational principle has been successfully applied to a wide range of problems in 
the theory of many-particle systems. The first application of Bogoliubov inequality to concrete 
many-particle problem was carried out in the work by I. A. Kvasnikov [157] on the application of 
a variational principle to the Ising model of ferromagnetism. 

Ising model [1711158] is defined by the following Hamiltonian T~L (i.e. energy functional of variables; 
in this case the ’’spins” Si = FI on the N sites of a regular lattice in a space of dimension d) 

N N 

U = - 72 Y1 ■ , U S ' S J - i,nil V .S', (129) 

i<j =1 i= 1 

Here, Jij play a role of ’’exchange constants”, H is a (normalized) magnetic held, involving an 
interpretation of the model to describe magnetic ordering in solids (M = Si is ’’magnetiza¬ 
tion”; hbHNL is the Zeeman energy, i.e. is the energy gained due to application of the held). 

The main task is to calculate statistical sum Z 


z = ^exp -(n/ey 

Si 

(130) 

Kvasnikov [157 considered the approximation of nearest neighbours, i.e. 
neighbours < i.j >. 

According to Bogoliubov variational principle one can write 

Jij = J for nearest 

F<F 0 + (H-no) 0 . 

(131) 

The upper bound for the free energy F sup is given by 


-^sup = 

(132) 
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where 


S = (Z 0 ) Tr( AH exp- 


(n 0 /e)i 

(133) 

[Bo m). 

(134) 


The parameters of partition, which were introduced into 7~Lq and A Ti, and, hence, into Zq and S, 
should be determined from the condition of the minimum of E sup . Thus we obtain 


Bp 


N 


a n 


N 


9 


= I j -bX 


2—1 


-x)^2 s i, 

(135) 

i=i 


N 


xi BijSiSj, 

(136) 




where B = H/9 , K = J/6 , x is some parameter. Then, according to relation (|132p . one finds 

-l 


Zinf 


(137) 
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2 coth hb{B - x) exp|/x s xtanh/r s (5 - x) + ^ ^ -fQj tanh 2 - x)}- 




Parameter x is determined by the equation 


tanh hb (x ~ B) = 


N 


E-K.J 


-/hBX, 


> 1 _ ^ XI K a + 


N 


E^ 


■(hsx ) 2 > 0 . 


(138) 


When the approximation of nearest neighbours is considered in the above equation the following 
substitution should be done 

N 

K ij = * KN , (139) 

where z is the number of nearest neighbours. 

Hence F sup is an approximate expression for the free energy and Z- m ^ is the approximate statistical 
sum of the model. It will be of instruction to compare these values with those which were calculated 
by other methods. To proceed, let us consider the regions of low and high temperatures. In the 
first case we will have that 9 <C zJ. The low-temperature approximation is expressed as a series 
expansion in terms of the small parameter exp(— K). The iterative solution of the Eg. (11381) will 
have the form 


I^bX = ~zK ^1 — 2 exp 2(—Kz — hbB) — 8 zK exp4 (—Kz — ^bB) + .. .^. (140) 

It is sufficient to confine oneself to the values of the order exp(— 2Kz). The result is 

(z-itfj = exp(Kz/2 + hbB)(i + exp 2(—Kz — B) + .. .^. (141) 

This result is in accordance with the other low-temperature expansions [5, 158 ) 

/ N z 

Z = exp 2(Kz/2 + hbB)N yl + N exp 2 {—Kz — fisB) -|—— exp 4 [—K(z — 1) — hbB] 

+{ N ^ 2 ^ ~ ex P 4 {-Kz - h b B) + ...)■ (142) 
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In the case of the high temperature, when 6 > zJ , the approximate solution of the 
have the form 

tanh bB 
zKj cosh 2 /jbB] 


PbX - -zK- 


Ea. (11381) will 
(143) 


Then after some transformations one can arrive to the expression (up to the terms it' 3 ): 


Z; n f ~ [2 cosh hbB] n ^1 + -Kztanh 2 hbB + (144) 

^K 2 zN[Az tanh 2 p bB + (Nz + 4 z) tanh 4 psBj'j . 

This expression is also in accordance with the known high-temperature expansions H3 IT58] for 
N » z. 

Let us consider now the expression for magnetization [[159] (the averaged magnetic moment) 


= I fllnZiM 
N dB ' 


(145) 


Using Eg. (11381) . we obtain 

m , /H m\ , . . 

-=tanh/x s (— +n —), (146) 

PbP v 0 0 / 

where p is the number of lattice sites per unit volume, m = Mp is the magnetization per unit vol¬ 
ume. This result coincides with the result of the phenomenological theory [5]. The corresponding 
basic values of the P. Weiss theory, the Curie point 6 o and Weiss parameter w have the form 
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y! Jij i 

*7O' 


w = 


N- 1 E Jr 


r-3 


V 2 B P 


Oq 

V 2 bP' 


(147) 


Hence, with the help of the Bogoliubov variational scheme it was possible to calculate the reason¬ 
able approximate expression for the statistical sum of the Ising model and describe the macroscopic 
properties of ferromagnetic systems in the wide interval of temperatures. It is thus seen that one 
may derive directly a consistent mean-field-type theory from a variational principle. 

Clearly Bogoliubov variational principle had a deep impact on the held of statistical mechanics of 
classical and quantum many-particle systems by making possible the analysis of complex statistical 
systems. Many interesting developments can be viewed from the point of a central theme, namely 
the Bogoliubov inequality, in particular in quantum theory of magnetism [5l I159L[160111611 162 . 163] 
and interacting many-body systems [164] 1165] 1166] 1167] 1168] 116911170111711 l9Tj . 

Radcliffe [l60j carried out a systematic investigation of the approximate free energies and Curie 
temperatures that can be obtained by using trial density matrices (which describe various possible 
decompositions of the ferromagnet into clusters) in a variational calculation of the free energy. 
Single-spin clusters lead to the molecular held model (as is well known) and two-spin clusters yield 
the Oguchi pair model m The relation of the constant-coupling method to these approximations 
was clarihed. A rigorous calculations using three-spin clusters were carried out. 

Rudoi [ 161] investigated the link between Bogoliubov statistical variational principle for free en¬ 
ergy, the method of partial diagram summation of the perturbation theory, and the Luttinger- 
Ward theorem. On the basis of Matsubara Green function method he solved the nonlinear integral 
Dyson equation by approximating the effective potential. As a result, a new implicit equation of 
magnetic state was obtained for the Ising model. 

Soldatov [172 ] generalized the Peierls-Bogoliubov inequality. A set of inequalities was derived in¬ 
stead, so that every subsequent inequality in this set approximates the quantity in question with 
better precision than the preceding one. These inequalities lead to a sequence of improving upper 
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bounds to the free energy of a quantum system if this system allows representation in terms of 
coherent states. It follows from the results obtained that nearly any upper bound to the ground 
state energy obtained by the conventional variational principle can be improved by means of the 
proposed method. 

F. Abubrig [173] studied the mixed spin-3/2 and spin-2 Ising ferrimagnetic system with different 
single-ion anisotropies in the absence of an external magnetic field within the mean-field the¬ 
ory based on Bogoliubov inequality for the Gibbs free energy. Second-order critical lines were 
obtained in the temperature-anisotropy plane. Tricritical line separating second-order and first- 
order lines was found. Finally, the existence and dependence of a compensation points on single-ion 
anisotropies was also investigated for the system. It was shown that this mixed-spin model exhibits 
one, two or three compensation temperature depending on the values of the anisotropies. 

10 The Variational Schemes and Bounds on Free Energy 

During last few decades, numerous variational schemes have become an increasingly popular work¬ 
able tool in quantum-mechanical many-particle theory mmmm- Bounds of free energy and 
canonical ensemble averages were of considerable interest as well. For many complex systems, 
such as Ising and Heisenberg ferromagnets or composite materials, methods of obtaining bounds 
are the practical useful tools which are both tractable and informative. A few illustrative topics 
will be of instruction to discuss in this context. 

MacDonald and Richardson [174]] used the density matrix of von Neumann to formulate an ex¬ 
act variational principle for quantum statistics which embodies the principle of maximization of 
entropy. In terms of the formalism of second quantization, authors wrote this variational princi¬ 
ple for fermions or bosons and then derived from it an approximate variational procedure which 
yields the particle states of a system of interacting bosons or fermions as well as the distribution 
of particles in these states. These equations, in authors’s opinion, yield the generalization of the 
Hartree-Fock equations for nonzero temperature and the corresponding extension to bosons. 

W. Schattke [175] found an upper bound for the free energy for superconducting system in mag¬ 
netic field. Starting from the BCS theory, the free energy was obtained by a combination of a 
variational method and perturbation theory. The variational equations obtained were non-local. 
The parameters of the perturbation calculation were the vector potential and the spatial varia¬ 
tions of the order parameter, which have to be small. Boundary conditions were set for the case 
of diffuse reflection and pair-breaking at the surface. As an example, the superconducting plate 
was discussed. 

Krinsky et. al. used [176] the variational principle to derive a new approximation to a ferromag- 
net in a magnetic field, below its critical temperature. They considered [176] a ferromagnet in an 
external magnetic field with T < Tq- Using a variational approximation based on the zero-field 
solution, the authors obtained an upper bound on the free energy, an approximate equation of 
state, and a lower bound on the magnetization, all having the correct critical indices. Explicit 
numerical calculations have been carried out for the two-dimensional Ising model, and it was 
found that the results obtained provide a good approximation to the results of series expansions 
throughout the region T < Tq. 

The Gibbs-Bogoliubov inequality [8] was used [177 ] to develop a first-order perturbation theory 
that provides an upper bound on the Helmholtz free energy per unit volume of a classical statisti¬ 
cal mechanical system in terms of the free energy and pair distribution function. Charged systems 
as well as a system of Lennard-Jones particles were discussed and detailed numerical estimates of 
the bounds were presented. 

S. Okubo and A. Isihara [178] derived important general inequalities for the derivatives of the par¬ 
tition function of a quantum system with respect to the parameters included in the Hamiltonian. 
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Applications of the inequalities were used to discuss relations for critical initial exponents, kinetic 
energy, susceptibility, electrical conductivity and so on. Existence of an inconsistency analogous 
to the Schwinger-term difficulty in the quantum held theory was pointed out. 

In their second paper [179] , S. Okubo and A. Isihara analyzed from a general point of view an 
inequality for convex functions in quantum-statistical mechanics. From an inequality for a convex 
function of two hermitian operators, the Peierls and Gibbs operators, coarse graining and other 
important inequalities were derived in a unified way. Various different forms of the basic inequality 
were given. They are found useful in discussing the entropy and other physical problems. Special 
accounts were given of functions such as exp(x) and xlogx. 

A variational method for many-body systems using a separation into a difference of Hamiltonians 
was presented by M. Hader and F. G. Mertens [180j . A particular ansatz for the wave function 
was considered which leads to an upper bound for the exact ground-state energy. This allowed a 
variation with respect to a separation parameter. The method was tested for a one-dimensional 
lattice with Morse interactions where the Toda subsystems can be solved by the Bethe ansatz. In 
two limiting cases results obtained were exact, otherwise they were in agreement with the quan¬ 
tum transfer integral method. 

Yeh [181, 118211183] proposed a derivation of a lower bound on the free energy; in addition he ana¬ 
lyzed the bounds of the average value of a function j l83j . He also established [ 181] a weaker form 
of Griffiths theorem for the ferromagnetic Heisenberg model. It was described as follows [182] : 
Free energy in the canonical ensemble was taken as 

F = — /3” 1 ln^(n| exp(— H/3)\n), (148) 

n 

where |n) is any complete set of orthonormal states. Bounds of F can be obtained from bounds 
of (n | exp(— Hj3)\n). As we seen, a very simple upper bound of F was given by Peierls [63]; one 
way to prove his theorem is by showing that 


(V ; |exp(— HP)\i!)) > exp(-/3(Y>|17|^)). 


(149) 


Yeh [182] derived a rather simple lower bound of F by similar method. He considered a a Hamil¬ 
tonian with a ground-state energy E$ = 0. He considered a real function f(E ) = exp(— E(3), 
(3 > 0. It was shown that for any normalized state | ip) a weaker but simpler upper bound for / 
may be written as 


exp 



(V’| exp (-17/3) |V>) > exp(-/3(V’|H#)), (150) 


where 


p = exp 




\H\ 




(151) 


Identifying /3 = (ksT) 1 and H as Hamiltonian, a lower bound of free energy was obtained from 
Eg. (11501) as 

F> -^ _i ln^exp(-/3(^|17|^)p), (152) 

where \i(j) is any complete orthonormal set of states. This is a general formula for a lower bound 
on the free energy. 

Upper and lower bounds of the canonical ensemble average of any operator A can be written down 
in terms of (ip n \H\ tp n ), where the <p n are eigenstates of A. Furthermore, bounds of thermodynamic 
derivatives can be obtained by noting that the bounds of 


d\f 

dj3 l 


( 153 ) 
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can be also derived [182] in similar manner. Here 


/ = (V’|exp(-H73)|V’) = ^2p n exp(-E n /3)- ^p n = l. (154) 

n n 

From Ea. (|150D . it is clear that all the bounds are more accurate at higher temperatures. These 
bounds have been useful in determining the properties of Heisenberg ferromagnets [181 ]. 

K. Symanzik [ 184 ] proved, refined, and generalized a lower bound given by Feynman for the quan¬ 
tum mechanical free energy of an oscillator. The method, application of a classical inequality to 
path integrals, also gives upper bounds for one-temperature Green functions. 

M. Heise and R. J. Jelitto [185 ] formulated the asymptotically exact variational approach to the 
strong coupling Hubbard model. They used a generalization of Bogoliubov variational principle, 
in order to develop a molecular field theory of the Hubbard model, which becomes asymptotically 
exact in the strong coupling limit. In other words, in their paper authors have started from a gen¬ 
eralized form of Bogoliubov variational theorem in order to set up a theory of the Hubbard model, 
which yields nontrivial results in the strong coupling regime and becomes asymptotically exact in 
the strong coupling limit. For this purpose the Hamiltonian was rotated by a unitary two-particle 
transformation, before the variational principle was applied. However, the real form of the gener¬ 
alized mean fields for the the Hubbard model in the strong coupling regime was not determined 
in complete form. This task was fulfilled by A. L. Kuzemsky in a series of papers [171 [35l [Ml [37]. 
K. Zeile [186 ] proposed a generalization of Feynman variational principle for real path integrals 
in a systematic way. He obtained an asymptotic series of lower bounds for the partition function. 
Author claimed that the method was tested on the anharmonic oscillator and showed excellent 
agreement with exact results. However, P. Dorre [187 ] et al. using the equivalence between Feyn¬ 
man and Bogoliubov variational principle, discussed [187] in the formalism of Hamiltonian quan¬ 
tum mechanics an improved upper bound for the free energy which has been given by Zeile [186] 
using path integral methods. It was shown that Zeile’s variational principle does not guarantee a 
thermodynamically consistent description. 

U. Brandt and J. Stolze formulated [ 188 ! a new hierarchy of upper and lower bounds on expecta¬ 
tion values. Upper and lower bounds were constructed for expectation values of functions of a real 
random variable with derivatives up to order [N + 1) which are alternately negative and positive 
over the whole range of interest. The bounds were given by quadrature formulas with weights 
and abscissas determined by the first (N + 1) moments of the underlying probability distribution. 
Application to a simple disordered phonon system yielded sharp bounds on the specific heat. 

K. Vlachos proposed a variational method that uses the frequency and the energy shift as 
variational parameters. The quantum-mechanical partition function was approximated by a for¬ 
mally simple expression, for a generalized anharmonic oscillator in one and many dimensions. The 
numerical calculations for a single quartic and two coupled quartic oscillators have led to nearly 
exact values for the free energy, the ground state, and the difference between the ground state 
and the first excited state. 

C. Predescu [ 189 ] presented a generalization of the Gibbs-Bogoliubov-Feynman inequality for 
spinless particles and then illustrated it for the simple model of a symmetric double-well quartic 
potential. The method gives a pointwise lower bound for the finite-temperature density matrix 
and it can be systematically improved by the Trotter composition rule. It was also shown to pro¬ 
duce ground state energies better than the ones given by the Rayleigh-Ritz principle as applied 
to the ground state eigenfunctions of the reference potentials. Based on this observation, it was 
conjectured that the local variational principle may performs better than the equivalent methods 
based on the centroid path idea and on the Gibbs-Bogoliubov-Feynman variational principle, es¬ 
pecially in the range of low temperatures. However clear evidence for such a statement was not 
given. 
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All these points of view acquire significance of the variational principles as a general method of 
solution for better insight into the complicated behavior of the many-particle systems. 

11 The Hartree-Fock-Bogoliubov Mean Fields 

The microscopic theory of superconductivity was created simultaneously by Bardeen, Cooper, and 
Schrieffer [190|. 1191] and Bogoliubov [19211193] 119411195] 11961119711198] . An important contribution 
to the theory of superconductivity were the works of Frohlich [199]. who put forward the idea of 
the importance of the electron-phonon interaction for the phenomenon of superconductivity, and 
the theory of Schafroth, Butler, and Blatt 2001 . who conjectured that superconductivity is due 
to Bose-Einstein condensation of correlated electron pairs. In their paper, Bardeen, Cooper, and 
Sehrieffer determined the ground state energy, and the spectrum of elementary excitations of their 
model [ 190L 191] . The BCS theory was constructed on the basis of a model Hamiltonian that takes 
into account only the interaction of electrons with opposite momenta and spins, whereas Bogoli¬ 
ubov theory was based on the Frohlich Hamiltonian [199] and used the method of compensation 
of dangerous diagrams [193] , N. N. Bogoliubov, V. V. Tolmachev and D. V. Shirkov [194] have 
generalized to Fermi systems the Bogoliubov method of canonical transformations proposed earlier 
in connection with a microscopic theory of superfluidity for Bose systems m- This approach 
has formed the basis of a new method for investigating the problem of superconductivity. Starting 
from Frohlich Hamiltonian, the energy of the superconducting ground state and the one-fermion 
and collective excitations corresponding to this state were obtained. It turns out that the final 
formulae for the ground state and one-fermion excitations obtained independently by Bardeen, 
Cooper and Schrieffer [190] were correct in the first approximation. The physical picture appears 
to be closer to the one proposed by Schafroth, Butler and Blatt. The effect on superconductivity 
of the Coulomb interaction between the electrons was analyzed in detail. A criterion for the su¬ 
perfluidity of a Fermi system with a four-line vertex Hamiltonian was established. 

Roughly speaking, to explain simply the theory of superconductivity it is possible to say that the 
Fermi sea is unstable against the formation of a bound Cooper pair when the net interaction is 
attractive; it is reasonable to expect that the pairs will be condense until an equilibrium point will 
be reached. The corresponding antisymmetric wave functions for many electrons was constructed 
in BCS model [ 201 . 202]. They noted also that their solution may be considered as an exact in 
the thermodynamic limit. 

The most clear and rigorous arguments in favor of the statement that the BCS model is an ex¬ 
actly solvable model of statistical physics were advanced in the papers of Bogoliubov, Zubarev, 
and Tserkovnikov [ 192 . 197 . 198 ], They showed that the free energy and the correlation functions 
of the BCS model and a model with a certain approximating quadratic Hamiltonian are indeed 
identical in the thermodynamic limit. In his theory [ 19211193111941 ! 19-3 . [196. 197 . 198] , Bogoliubov 
gave a rigorous proof that at vanishing temperature the correlation functions and mean values of 
the energy of the BCS model and the Bogolyubov-Zubarev-Tserkovnikov model are equal in the 
thermodynamic limit. Moreover, Bogoliubov constructed a complete theory of superconductivity 
on the basis of a model of interacting electrons and phonons [192 . 1193 1 11941 195 . 11961119711198] . 
Generalizing his method of canonical transformations [151 l2D3l 204] to Fermi systems and ad¬ 
vancing the principle of compensation of dangerous graphs [194] , he determined the ground state 
consisting of paired electrons with opposite moments and spins, its energy, and the energy of 
elementary excitations. It was shown also that the phenomenon of superconductivity consists 
in the pairing of electrons and a phase transition from a normal state with free electrons to a 
superconducting state with pair condensate. 
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The pairing Hamiltonian has the form 

H- pM = E E{k)a\ a a ka + ^V(k, p)a\^_ ki a- pi a p ^, (155) 

ka kp 

where p is the chemical potential and A/" is the number of particles. 

The essential step which was made by Bogoliubov was connected with introducing the anomalous 
averages or the generalized mean fields F p = (a- p ±a p f). It is reasonably to suppose that because 
of the large number of particles involved, the fluctuations of a- p ^a p ^ about these expectations 
values F p must be small. Hence, it is thus possible to express such products of operators in the 
form 

= Fp F (®— F p )- (156) 

It is reasonably to suppose that one may neglects by the quantities which are bilinear in the pre¬ 
sumably small fluctuation term in brackets. This way leads to the Bogoliubov model Hamiltonian 
of the form 


T~Lmod phf — 

E ( k ) a la a k<r + Y V ( k ’P) { a U a -ki F p + F^a_ pi a pt - FfF^j . (157) 

ka kp 

Here the F k should be determined self-consistently [19211193111941119511196111971 [l98j . 

Thus Bogoliubov created a rigorous theory of superfluidity [112] and superconductivity [198] within 
the unified scheme [115] of the non-zero anomalous averages or the generalized mean fields, and 
showed that at the physical basis of these two fundamental phenomena of nature lies the process 
of condensation of Bose particles [1T6] and, respectively, pairs of fermions. 

Indeed, N. N. Bogoliubov, D. N. Zubarev and Yu. A. Tserkovnikov [1921 [197] , have shown on 
the basis of the model Hamiltonian of BCS-Bogoliubov, that the thermodynamic functions of a 
superconducting system, which were obtained by a variation method in BCS, are asymptotically 
exact for V —>• oo ,N/V = const (V is the volume of the system, and N the number of parti¬ 
cles). This conclusion was based on the fact that each term of the perturbation-theory series, 
by means of which the correction to that solution was calculated, is asymptotically small for 
V —>• oo. In addition, it was shown that it is possible to satisfy with asymptotic exactness the 
entire chain of equations for Green functions constructed on the basis of the model Hamiltonian of 
BCS-Bogoliubov. Thus the asymptotic exactness of the known solution for the superconducting 
state was proved without the use of perturbation theory. It was shown also that the trivial solu¬ 
tion that corresponds to the normal state should be rejected at temperatures below the critical 
temperature. In other words, starting with the reduced Hamiltonian of superconductivity theory, 
Bogoliubov, Zubarev, and Tserkovnikov .192.1197] proved the possibility of exact calculation of 
the free energy per unit volume. 

Somewhat later, on the basis of the BCS theory, a similar investigation was made by other au¬ 
thors [2051 12061 12071 1208] . B. Muhlschlegel [205] studied an asymptotic expansion of the BCS 
partition function by means of the functional method. The canonical operator exp [—/3(H — pN)] 
associated with the BCS model Hamiltonian of superconductivity was represented as a functional 
integral by the use of Feynman’s ordering parameter. General properties of the partition function 
in this representation were investigated. Taking the inverse volume of the system as an expansion 
parameter, it was possible to calculate the thermodynamic potential including terms independent 
of the volume. Muhlschlegel’s theory yielded an additional evidence that the BCS variational 
value is asymptotically exact. The behavior of the canonical operator for large volume was de¬ 
scribed and related to the state of free quasiparticles. A study of the terms of the thermodynamic 
potential which were of smaller order in the volume in the low-temperature limit, showed that the 
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ground state energy is nondegenerate and belongs to a number eigenstate. 

W. Thirring and A. Wehrl [209] investigated in which sense the Bogoliubov-Haag treatment of 
the BCS-Bogoliubov model gives the correct solution in the limit of infinite volume. They found 
that in a certain subspace of the infinite tensor product space the field operators show the correct 
time behaviour in the sense of strong convergence. Thus a solution of the superconducting type 
with a gap in the spectrum of elementary excitations really can exists for the model Hamiltonian 
of BCS-Bogoliubov. 

In general, the problem of explaining the phenomenon of superconductivity required the solution 
of the very difficult mathematical problems associated with the foundation of applied approxima¬ 
tions PlUg. In connection with this, Bogoliubov investigated |192l 11931119411195111961119711198] 
the reduced Hamiltonian, in which the interaction of single electrons is studied, and carried out 
for it a complete mathematical investigation for zero temperature. In this connection he laid the 
bases of a new powerful method of the approximating Hamiltonian , which allows linearization 
of nonlinear quantum equations of motion, and reduction of all nonlinearity to self-consistent 
equations for the ordinary functions into which the defined operator expressions translate. This 
method was extended later to nonzero temperatures and a wide class of systems, and became one 
of most powerful methods of solving nonlinear equations for quantum fields mm- 
D. Ya. Petrina contributed much to the further clarification of many complicated aspects of the 
BCS-Bogoliubov theory. He performed a close and subtle analysis [ T5l [TT5112101 [ZTI , [2T2, 12T3] of 
the BCS-Bogoliubov model and various related mathematical problems. 

In his paper [210] ” Hamiltonians of quantum statistics and the model Hamiltonian of the theory 
of superconductivity” an investigation was made of the general Hamiltonian of quantum statistics 
and the model Hamiltonian of the theory of superconductivity in an infinite volume. The Hamilto¬ 
nians were given a rigorous mathematical definition as operators in a Hilbert space of sequences of 
translation-invariant functions. It was established that the general Hamiltonian is not symmetric 
but possesses a real spectrum; the model Hamiltonian is symmetric and its spectrum has a gap 
between the energy of the ground state and the excited states. 

In the following paper [[21 ij . the model Hamiltonian of the theory of superconductivity was in¬ 
vestigated for an infinite volume and a complete study was made of its spectrum. The grand 
partition function was determined and the equation of state was found. In addition, the existence 
of a phase transition from the normal to the superconducting state was proved. It was shown 
that in the limit V —> oo the chain of equations for the Green functions of the model Hamiltonian 
has two solutions, namely the free Green function and the Green function of the approximating 
Hamiltonian. 

In his paper [212] D.Petrina has shown that the Bogolyubov result that the average energies (per 
unit volume) of the ground states for the BCS-Bogoliubov Hamiltonian and the approximating 
Hamiltonian asymptotically coincide in the thermodynamic limit is also valid for all excited states. 
He also established that, in the thermodynamic limit, the BCS-Bogoliubov Hamiltonian and the 
approximating Hamiltonian asymptotically coincide as quadratic forms. 

D.Petrina [213] considered also the BCS Hamiltonian with sources, as it was proposed by Bo¬ 
goliubov and Bogoliubov, Jr. It was proved that the eigenvectors and eigenvalues of the BCS- 
Bogoliubov Hamiltonian with sources can be exactly determined in the thermodynamic limit. 
Earlier, Bogoliubov proved that the energies per volume of the BCS-Bogoliubov Hamiltonian 
with sources and the approximating Hamiltonian coincide in the thermodynamic limit. These 
results clarified substantially the microscopic theory of superconductivity and provided a deeper 
mathematical foundation to it. 

Raggio and Werner [ 169] have shown the existence of the limiting free-energy density of inhomoge¬ 
neous (site-dependent coupling) mean-field models in the thermodynamic limit [45], and derived 
a variational formula for this quantity. The formula requires the minimization of an energy term 
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plus an entropy term as a functional depending on a function with values in the one-particle state 
space. The minimizing functions describe the pure phases of the system, and all cluster points 
of the sequence of finite volume equilibrium states have unique integral decomposition into pure 
phases. Some applications were considered; they include the full BCS-model, and random mean- 
field models. 

A detailed and careful mathematical analysis of certain aspects of the BCS-Bogoliubov theory was 
carried out by S. Watanabe [2141121511216112171121811219112201122111222] , mainly in the context of 
the solutions to the BCS-Bogoliubov gap equation for superconductivity. 

BCS-Bogoliubov theory correctly yields an energy gap [223, 224 ], The determination of this im¬ 
portant energy gap is by solving a nonlinear singular integral equation. An investigation of the 
solutions to the BCS-Bogoliubov gap equation for superconductivity was carried out by S. Watan¬ 
abe [2141121511216112171121811219112201122111222] . In his works the BCS-Bogoliubov equations were 
studied in full generality. Watanabe investigated the gap equation in the BCS-Bogoliubov theory 
of superconductivity, where the gap function is a function of the temperature T only. It was shown 
that the squared gap function is of class C 2 on the closed interval [0, Tc]. Here, Tc stands for 
the transition temperature. Furthermore, it was shown that the gap function is monotonically 
decreasing on [0, Tc] and the behavior of the gap function at T = Tc was obtained and some 
more properties of the gap function were pointed out. 

On the basis of his study Watanabe then gave, by examining the thermodynamical potential, a 
mathematical proof that the transition to a superconducting state is a second-order phase tran¬ 
sition. Furthermore, he obtained a new and more precise form of the gap in the specific heat 
at constant volume from a mathematical point of view. It was shown also that the solution to 
the BCS-Bogoliubov gap equation for superconductivity is continuous with respect to both the 
temperature and the energy under the restriction that the temperature is very small. Without 
this restriction, the solution is continuous with respect to both the temperature and the energy, 
and, moreover, the solution is Lipschitz continuous and monotonically decreasing with respect to 
the temperature. 

D. M. van der Walt, R. M. Quick and M. de Llano [225ll226 j have obtained analytic expressions for 
the BCS-Bogoliubov gap of a many-electron system within the BCS model interaction in one, two, 
and three dimensions in the weak coupling limit, but for arbitrary interaction width is = HD/Ep, 
0 < is < oo. Here HD is the maximum energy of a force-mediating boson, and Ep is the Fermi 
energy (which is fixed by the electronic density). The results obtained addressed both phononic 
(is -C 1) as well as nonphononic (e.g., exciton, magnon, plasmon, etc.) pairing mechanisms where 
the mediating boson energies are not small compared with Ep, provided weak electron-boson 
coupling prevails. The essential singularity in coupling, sometimes erroneously attributed to the 
two-dimensional character of the BCS model interaction with (is <C 1), was shown to appear in 
one, two, and three dimensions before the limit is —> 0 is taken. 

B. McLeod and Yisong Yang [227] studied the uniqueness and approximation of a positive solu¬ 
tion of the BCS-Bogolyubov gap equation at finite temperatures. When the kernel was positive 
representing a phonon-dominant phase in a superconductor, the existence and uniqueness of a 
gap solution was established in a class which contains solutions obtainable from bounded domain 
approximations. The critical temperatures that characterize superconducting-normal phase tran¬ 
sitions realized by bounded domain approximations and full space solutions were also investigated. 
It was shown under some sufficient conditions that these temperatures are identical. In this case 
the uniqueness of a full space solution follows directly. Authors [ 227] also presented some exam¬ 
ples for the nonuniqueness of solutions. The case of a kernel function with varying signs was also 
considered. It was shown that, at low temperatures, there exist nonzero gap solutions indicating 
a superconducting phase, while at high temperatures, the only solution is the zero solution, rep¬ 
resenting the dominance of the normal phase, which establishes again the existence of a transition 
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temperature. 

In a series of papers [228 . 22!) 230 . M. Combescot, W.V. Pogosov and O. Betbeder-Matibet stud¬ 
ied various aspects of the BCS ansatz for superconductivity [ 190] in the light of the Bogoliubov 
approach. 

In paper |228) they extended the one-pair Cooper configuration towards BCS-Bogoliubov model of 
superconductivity by adding one-by-one electron pairs to an energy layer, where a small attraction 
acts. To do it, they solved Richardson’s equations analytically in the dilute limit of pairs on the 
one-Cooper pair scale. It was found, through only keeping the first order term in this expansion, 
that the N correlated pair energy reads as the energy of N isolated pairs within a N(N — 1) cor¬ 
rection induced by the Pauli exclusion principle which tends to decrease the average pair binding 
energy when the pair number increases. Quite remarkably, extension of this first-order result to 
the dense regime gives the BCS-Bogoliubov condensation energy exactly. These facts may lead 
ones to a different interpretation of the BCS-Bogoliubov condensation energy with a pair number 
equal to the number of pairs feeling the potential and an average pair binding energy reduced by 
Pauli blocking to half the single Cooper pair energy - instead of the more standard but far larger 
superconducting. 

In the next paper [229] the usual formulation of the BCS-Bogoliubov ansatz for superconductiv¬ 
ity in the grand canonical ensemble makes the handling of the Pauli exclusion principle between 
paired electrons straightforward. It however masks that many-body effects between Cooper pairs 
interacting through the reduced BCS-Bogoliubov potential are entirely controlled by this exclu¬ 
sion. To show it up, one has to work in the canonical ensemble. The proper handling of Pauli 
blocking between a fixed number of composite bosons is however known to be quite difficult. To 
do it, author have developed a commutator formalism for Cooper pair condensate, along the line 
they used for excitons. Authors [229 ] then rederived, within the IV-pair subspace, a few results 
of BCS-Bogoliubov theory of superconductivity obtained in the grand canonical ensemble, to evi¬ 
dence their Pauli blocking origin. They ended by reconsidering what should be called Cooper pair 
wave function and concluded differently from usual understanding. 

In their third paper M. Combescot, W.V. Pogosov and O. Betbeder-Matibet [230] showed that the 
Bogoliubov approach to superconductivity provides a strong mathematical support to the wave 
function ansatz proposed by Bardeen, Cooper and Schrieffer [190 ]. However there are some subtle 
differences in the both the approaches. Indeed, the BCS ansatz - with all pairs condensed into the 
same state - corresponds to the ground state of the Bogoliubov Hamiltonian. From the other hand, 
this Hamiltonian only is part of the BCS Hamiltonian. As a result, the BCS ansatz definitely 
differs from the BCS Hamiltonian ground state. This can be directly shown either through a per¬ 
turbative approach starting from the Bogoliubov Hamiltonian, or better by analytically solving 
the BCS Schrodinger equation along Richardson-Gaudin exact procedure. Still, the BCS ansatz 
leads not only to the correct extensive part of the ground state energy for an arbitrary number 
of pairs in the energy layer where the potential acts - as recently obtained by solving Richardson- 
Gaudin equations analytically - but also to a few other physical quantities such as the electron 
distribution, as it was shown by authors. The paper [230] also considered arbitrary filling of the 
potential layer and evidences the existence of a super dilute and a super dense regime of pairs, 
with a gap different from the usual gap. These regimes constitute the lower and upper limits of 
density-induced BEC-BCS cross-over in Cooper pair systems. It should be noted, however, that 
this theory needs an additional careful examination. 

In 1958 N. N. Bogoliubov [231 ] proposed a new variational principle in the many-particle problem. 
This variational principle is the generalization of the Hartree-Fock variational principle [5][10]. It 
is well known [2321 [233] that the Hartree-Fock approximation is a variational method that pro¬ 
vides the wave function of a many-body system assumed to be in the form of a Slater determinant 
for fermions and of a product wave function for bosons. It treats correctly the statistics of the 
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many-body system, antisymmetry for fermions and symmetry for bosons under the exchange of 
particles. The variational parameters of the method are the single-particle wave functions com¬ 
posing the many-body wave function. 

Bogoliubov [231 ] considered a model dynamical Fermi system describing by the Hamiltonian with 
two-body forces. The Hamiltonian of a nonrelativistic system of identical fermions interacting by 
two-body interactions was 

H = ^ ( E(k ) - Ep'j a\ a a ka + 7^7 22 J(k, k'\aia 2 (r' 2 a' l )a\ (T a\ a akaaka- (158) 

ka /c,/c',cr 

The at and aka are single-particle creation and annihilation operators satisfying the usual anti¬ 
commutation relations, Ef is the Fermi energy level and V is the volume of the system. 

The Hamiltonian under consideration is a model Hamiltonian; it takes into account the pair in¬ 
teraction of the particles with opposite momentum only. It can be rewritten in the following 
form [231] : 

H = 22 ( E ( k ) ~ E r) a \s a qs + 7 ^ 7 / (9>9 , l a i> a 2,44)<4u 0 k'V4 a <?'sT ( 159 ) 

i s s 

Here q describes the pair of momentum (k, — k); hence q and — q describe the same pair. Index 
s = (cTju), where v = ±1 is an additional index [231 ] permitting to classify k as (q,v). N. N. 
Bogoliubov [231 ] shown that the ground state of the system can be found asymptotically exactly 
for the limit V —>• 00 by following to approach of the paper [ 192 ]. 

This approach found numerous applications in the many-body nuclear theory [232, 233 . 234 . 235 . 
1236 ; (23711238112391 [240] . The properties of all existing and theoretically predicted nuclei can be 
calculated based on various nuclear many-body theoretical frameworks. The classification of nu¬ 
clear many-body methods can be also done from the point of view of the pair nuclear interaction, 
from which the many-body Hamiltonian is constructed. An important goal of nuclear structure 
theory is to develop the computational tools for a systematic description of nuclei across the 
chart of the nuclides. Nuclei come in a large variety of combinations of protons and neutrons 
(< 300). Understanding the structure of the nucleus is a major challenge. To study some collec¬ 
tive phenomena in nuclear physics, we have to understand the pairing correlation due to residual 
short-range correlations among the nucleons in the nucleus. This has usually been calculated 
by using the BCS theory or the Hartree-Fock-Bogoliubov theory. The Hartree-Fock-Bogoliubov 
theory is suited well for describing the level densities in nuclei. [23711239] , The theory of level 
densities reminds in certain sense the ordinary thermodynamics. The simplest level density of 
nucleons calculations were based usually on a model Hamiltonian which included a simple version 
of the pairing interaction (between nucleons in states differing only by the sign of the magnetic 
quantum number). 

J. A. Sheikh and P. Ring [236] derived the symmetry-projected Hartree-Fock-Bogoliubov equa¬ 
tions using the variational ansatz for the generalized one-body density-matrix in the Valatin form. 
It was shown that the projected-energy functional can be completely expressed in terms of the 
Hartree-Fock-Bogoliubov density-matrix and the pairing-tensor. The variation of this projected- 
energy was shown to result in Hartree-Fock-Bogoliubov equations with modified expressions for 
the pairing-potential and the Hartree-Fock field. The expressions for these quantities were ex¬ 
plicitly derived for the case of particle number-projection. The numerical applicability of this 
projection method was studied in an exactly soluble model of a deformed single-j shell. 

A. N. Behkami and Z. Kargar [237 ] have determined the nuclear level densities and thermodynamic 
functions for light A nuclei, from a microscopic theory, which included nuclear pairing interaction. 
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Nuclear level densities have also been obtained using Bethe formula as well as constant tempera¬ 
ture formula. Level densities extracted from the theories were compared with their corresponding 
experimental values. It was found that the nuclear level densities deduced by considering various 
statistical theories are comparable; however, the Fermi-gas formula [241 ] becomes inadequate at 
higher excitation energies. This conclusion, which has also been arrived at by other investigations, 
revealed that a realistic treatment of the statistical nuclear properties requires the introduction 
of residual interaction. The effects of the pairing interaction and deformation on nuclear state 
densities were illustrated and discussed. 

L. M. Robledo and G. F. Bertsch [ 238] have presented a computer code for solving the equa¬ 
tions of the Hartree-Fock-Bogoliubov theory by the gradient method, motivated by the need for 
efficient and robust codes to calculate the configurations required by extensions of the Hartree- 
Fock-Bogoliubov theory, such as the generator coordinate method. The code was organized with a 
separation between the parts that are specific to the details of the Hamiltonian and the parts that 
are generic to the gradient method. This permitted total flexibility in choosing the symmetries 
to be imposed on the Hartree-Fock-Bogoliubov solutions. The code solves for both even and odd 
particle-number ground states, with the choice determined by the input data stream. 

M. Lewin and S. Paul [240 ] shown that the best method for describing attractive quantum systems 
is the Hartree-Fock-Bogoliubov theory. This approach deals with a nonlinear model which allows 
for the description of pairing effects, the main explanation for the superconductivity of certain 
materials at very low temperature. Their paper is a detailed study of Hartree-Fock-Bogoliubov 
theory from the point of view of numerical analysis. M. Lewin and S. Paul started by discussing 
its proper discretization and then analyzed the convergence of the simple fixed point (Roothaan) 
algorithm. Following works for electrons in atoms and molecules, they shown that this algorithm 
either converges to a solution of the equation, or oscillates between two states, none of them being 
solution to the Hartree-Fock-Bogoliubov equations. They also adapted the Optimal Damping Al¬ 
gorithm to the Hartree-Fock-Bogoliubov setting and also analyzed it. The last part of the paper 
was devoted to numerical experiments. Authors considered a purely gravitational system and 
numerically discovered that pairing always occurs. They then examined a simplified model for 
nucleons, with an effective interaction similar to what is often used in nuclear physics. In both 
cases M. Lewin and S. Paul [240] discussed the importance of using a damping algorithm. 

Many other applications of the Hartree-Fock-Bogoliubov theory to various many-particle systems 
were discussed in Refs. [24211243[ [244112451 [2461 Generalization of Lieb variational principle [ 166] 
to Bogoliubov-Hartree-Fock theory was considered recently by V. Bach et al. [ 167 ] In its original 
formulation, Lieb variational principle holds for fermion systems with purely repulsive pair interac¬ 
tions. As a generalization authors proved for both fermion and boson systems with semi-bounded 
Hamiltonian that the infimum of the energy over quasifree states coincides with the infimum over 
pure quasifree states. In particular, the Hamiltonian was not assumed to preserve the number of 
particles. 

It is instructive to remind that in mathematics, the infimum (abbreviated inf; plural infima) of 
a subset S of a partially ordered set T is the greatest element of T that is less than or equal to 
all elements of S. Consequently the term greatest lower bound is also commonly used. Infima 
of real numbers are a common special case that is especially important in analysis. However, 
the general definition remains valid in the more abstract setting of order theory where arbitrary 
partially ordered sets are considered. 

To shed light on the relation between authors’ result and the usual formulation of Lieb variational 
principle in terms of one-particle density matrices, it was also included a characterization of pure 
quasifree states by means of their generalized one-particle density matrices. 
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12 Method of an Approximating Hamiltonian 


It is worth noting that a complementary method, which was called by the method of an approxi¬ 
mating Hamiltonian , was formulated [3], 0] '24711248 - ; 239] for treating model systems of statistical 
mechanics. The essence of the method consists in replacement of the initial model Hamiltonian 
H , which is not amenable to exact solution, by a suitable approximating (or trial) Hamiltonian 
H 3 P pr . The next step consists of proving their thermodynamical equivalence, i.e., proving that 
the thermodynamic potentials and the mean values calculated on the basis of H and H appr are 
asymptotically equal in the thermodynamic limit [45| N, V —>• oo, N/V = const. 

When investigating the phenomenon of superconductivity, Bogolyubov suggested the method of 
approximating Hamiltonian and justified it for the case of temperatures close to zero. By employ¬ 
ing this method, Bogolyubov rigorously solved the BCS model of superconductivity at temperature 
zero. This model was defined by the Hamiltonian of interacting electrons with opposite momenta 
and spins. 

To explain the superconductivity phenomenon, it was necessary to solve very difficult mathemat¬ 
ical problems connected with the justification of approximations employed. In this connection, 
Bogoliubov considered the reduced Hamiltonian in which only the interaction of electrons was 
taken into account. He gave a complete mathematical investigation of this Hamiltonian at tem¬ 
perature zero. Moreover, he laid the foundation of a new powerful method of approximating 
Hamiltonian which allows one to linearize nonlinear quantum equations of motion so that the 
nonlinearity is preserved only in self-consistent equations for ordinary functions that are obtained 
from certain operator expressions. This method was then extended to the case of nonzero tem¬ 
peratures and applied to a broad class of systems. Later, this approach became one of the most 
effective methods for solving nonlinear equations for quantum fields. 

The method of approximating Hamiltonian is based on the proof of the thermodynamic equivalence 
of the model under consideration and approximating Hamiltonian. Thermodynamic equivalence 
means here the coincidence of specific free energies and Green functions for model and approxi¬ 
mating Hamiltonian in the thermodynamic limit [35] when V and N tends to oo, N/V = const. 
It was shown above that in many cases it may be assumed that the effective Hamiltonian H for the 
system of particles may be written as the sum of the Hamiltonian of the reference system i7 appr , 
plus the rest of the effective Hamiltonian H = 77 appr + AH. Then the Bogoliubov inequality 
states that the Helmholtz free energy F of the system is given by 

F < T appr + (H- LT appr ) appr , (160) 

where I pappr notes the free energy of the reference system and the brackets a canonical ensemble 
average over the reference system. 

N.N. Bogoliubov Jr. elaborated a new method (247 . 1248 . 249 . 250 . 251 ] of finding exact solutions 
for a broad class of model systems in quantum statistical mechanics - the method of approximating 
Hamiltonian. As it was mentioned above this method appeared in the theory of superconductiv¬ 
ity [1971IT98[ . 

N.N. Bogoliubov Jr. investigated some dynamical models (247 ] generalizing those of the BCS 
type. A complete proof was presented that the well-known approximation procedure leads to an 
asymptotically exact expression for the free energy, when the usual limiting process of statistical 
mechanics is performed. Some special examples were considered. 

A detailed analysis of Bogoliubov approach to investigations of (Hartree-Fock-Bogoliubov) mean- 
field type approximations for models with a four-fermion interaction was given in the papers (250 . 
1251] . An exactly solvable model with paired four-fermion interaction that is of interest in the 
theory of superconductivity was considered. Using the method of approximating Hamiltonian, 
it was shown that it is possible to construct an asymptotically exact solution for this model. 
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In addition a theorem was proved that allows us to compute, with asymptotic accuracy in the 
thermodynamic limit, the density of the free energy under sufficiently general conditions imposed 
on the parameters of the model system. An approximate method for investigating models with 
four-fermion interaction of general form was presented. The method was based on the idea of 
constructing an approximating Hamiltonian and it allows one to study the dynamical properties 
of these models. The method combines the standard approach to the method of the approximat¬ 
ing Hamiltonian for the investigation of models with separable interaction and the Hartree-Fock 
scheme of approximate computations based on the concept of self-consistency. To illustrate the 
efficiency of the approach presented, the BCS model that plays an important role in the theory of 
superconductivity was considered in detail. Thus, the effective and workable approach was formu¬ 
lated which allows one to investigate dynamical and thermodynamical properties of models with 
four-fermion interaction of general type. The approach combines ideas of the standard Bogoliubov 
approximating Hamiltonian method for the models with separable interaction with the method of 
Hartree-Fock approximation based on the ideas of self-consistency. 

A. P. Bakulev, N. N. Bogoliubov, Jr., and A. M. Kurbatov [252 ] discussed thoroughly the prin¬ 
ciple of thermodynamic equivalence in statistical mechanics in the approach of the method of 
approximating Hamiltonian. They discussed the main ideas that lie at the foundations of the 
approximating Hamiltonian method in statistical mechanics. The principal constraints for model 
Hamiltonian to be investigated by approximating Hamiltonian method were considered along with 
the main results obtainable by this method. It was shown how it is possible to enlarge the class 
of model Hamiltonians solvable by approximating Hamiltonian method with the help of an ex¬ 
ample of the BCS-type model. Additional rigorous studies of the theory of superconductivity 
with Coulomb-like repulsion was carried out by A. P. Bakulev [253] . The traditional method of 
the approximating hamiltonian was applied for the investigation of a model of a superconductor 
with interaction of the BCS type and Coulomb-like repulsion, the latter being described by un¬ 
bounded operators. It was shown that the traditional method can be generalized in such a way 
that for the model under consideration one can prove the asymptotic (in the thermodynamic limit 
V —>• oo,N —» oo, N/V = const) coincidence not only of the free energies (per unit volume) but 
also of the correlation functions of the model and approximating Hamiltonian. 

13 Conclusions 

The aim of the present overview was to justify a statement that in many cases the methods of 
quantum statistical mechanics, many of which were formulated and developed by N.N. Bogoli¬ 
ubov mm, allow one to develop efficient approaches for solution of complicated problems of 
the many-particle interacting systems. 

In the present survey we discussed tersely the Bogoliubov variational principle. It was shown in 
the preceding sections that this principle provides an extremely valuable treatment of mean-field 
methods and their application to the problems in statistical mechanics and many-particle physics 
of interacting systems. With its remarkable workability the Bogoliubov variational principle found 
many applications as an effective method not only in condensed matter physics but also in many 
other areas of physics (see, e.g. Ref. [254]). It is also hoped that this work will lead to greater 
insight into the application of variational principles to various many-particle problems. 

There is another aspect of the problem under consideration. It is of great importance to determine 
correctly the mean-field contribution when one describes the interacting many-particle systems 
by the equations-of-motion method 13 DEI. It was mentioned briefly that the method of two-time 
temperature Green functions 13 E| allows one to investigate efficiently the quasiparticle rnany- 
body dynamics generated by the main model Hamiltonians from the quantum solid state theory 
and the quantum theory of magnetism. The method of quasiaverages allows one to take a deeper 
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look at the problems of spontaneous symmetry breaking, as well as at the problems of symmetry 
and dissymmetry in the physics of condensed matter [sum ns]. 

Summarizing the basic results obtained by N. N. Bogoliubov by inventing the variational prin¬ 
ciples, method of quasiaverages and results in the area of creation of asymptotic methods of 
statistical mechanics, one must especially emphasize that thanks to their deep theoretical con¬ 
tent and practical direction, these methods have obtained wide renown everywhere. They have 
enriched many-particle physics and statistical mechanics with new achievements in the area of 
mathematical physics as well as in the areas of concrete applications to physics, e.g. theories of 
superfluidity and superconductivity. 

In the papers [I7i SSISSIST], we have formulated the self-consistent theory of the correlation effects 
for many-particle interacting systems using the ideas of quantum field theory for interacting elec¬ 
tron and spin systems on a lattice. The workable and self-consistent irreducible Green functions 
approach to the decoupling problem for the equation-of-motion method for double-time temper¬ 
ature Green functions has been presented. The main achievement of this formulation was the 
derivation of the Dyson equation for double-time retarded Green functions instead of causal ones. 
That formulation permitted to unify convenient analytical properties of retarded and advanced 
Green functions and the formal solution of the Dyson equation, that, in spite of the required ap¬ 
proximations for the self-energy, provides the correct functional structure of single-particle Green 
function. The main advantage of the mathematical formalism was brought out by showing how 
elastic scattering corrections (generalized mean fields) and inelastic scattering effects (damping 
and finite lifetimes) could be self-consistently incorporated in a general and compact manner. We 
have presented there the novel method of calculation of quasi-particle spectra for basic spin lat¬ 
tice models, as the most representative examples. Using the irreducible Green functions method, 
we were able to obtain a closed self-consistent set of equations determining the electron Green 
function and self-energy. For the Hubbard and Anderson models, these equations gave a general 
microscopic description of correlation effects both for the weak and strong Coulomb correlation, 
and, thus, determined the interpolation solutions of the models. Moreover, this approach gave 
the workable scheme for the definition of relevant generalized mean fields written in terms of ap¬ 
propriate correlators. 

We hope that these methods of statistical mechanics have been explained with sufficient details 
to bring out their scope and power, since we believe that those techniques will have application 
to a variety of many-body systems with complicated spectra and strong interaction. 

These applications have illustrated some of subtle details of the irreducible Green functions ap¬ 
proach and exhibited their physical significance in a representative form. 

As it was seen, these treatments has advantages in comparison with the standard methods of 
decoupling of higher-order Green functions within the equation-of-motion approach. 

The main advantage of the whole method is the possibility of a self-consistent description of quasi¬ 
particle spectra and their damping in a unified and coherent fashion. 

The most important conclusion to be drawn from the present consideration is that the generalized 
mean fields for the case of strong Coulomb interaction in the Hubbard model has quite a nontrivial 
structure and cannot be reduced to the mean-density functional. 

Recently the problem of the advanced mean field methods in complex systems |255 ] has attracted 
big attention. Our consideration reveals the fundamental importance of the adequate definition of 
generalized mean fields at finite temperatures, that results in a deeper insight into the nature of 
quasiparticle states of the correlated lattice fermions and spins and other interacting many-particle 
systems. 
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